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1.0 Introduction 

Since the first extravehicular activities (EVAs) were per- 
formed by Alexei Leonov and Ed White in 1965, the 
capabilities of astronauts to do useful work outside of 
their spacecraft have steadily increased. Likewise, our 
understanding of EVA astronauts’ capabilities and limi- 
tations have also progressed through inflight experience, 
experimentation in neutral buoyancy facilities and para- 
bolic flight, and engineering tests of space suits and 
EVA tools. The purpose of this research effort is to fur- 
ther advance the current understanding of astronauts’ 
capabilities and limitations in space-suited EVA. 

The most important aspect of an EVA astronaut’s capa- 
bilities is the ability to move his or her body while wear- 
ing the space suit. In every EVA scenario, astronauts 
physically manipulate objects to accomplish tasks. Two 
factors make these physical interactions strikingly dif- 
ferent from those performed on the ground. First, the 
microgravity environment requires an astronaut to 
restrain his or her body in order to exert forces and 
moments on another object, and second, space suits con- 
strain astronauts’ body motions in significant and com- 
plicated ways. 

Historically, feasible limits on planned EVA activities 
have been determined based on ground experimentation 
and in-flight experience. Data have been obtained exper- 
imentally on operational performance metrics such as 
joint ranges of motion, the volume in space that a space- 
suited astronaut can reach, known as the reach envelope, 
the subset of the reach envelope in which a space-suited 
astronaut can comfortably work, known as the work 
envelope, and the strength of a suited astronaut. These 
performance measures are related to each other through 
the constitutive and compatibility relations that govern 
the space suit’s behavior, but isolated experimental data 
does not allow one performance measure to be predicted 
from another or the same perfonnance measure to be 
generalized to a different person or situation. 
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The primary aim of this research effort was to advance 
the current understanding of astronauts’ capabilities and 
limitations in space-suited EVA by developing models 
of the constitutive and compatibility relations of a space 
suit, based on experimental data gained from human test 
subjects as well as a 12 degree-of-ffeedom human-sized 
robot, and utilizing these fundamental relations to esti- 
mate a human factors performance metric for space 
suited EVA work. The three specific objectives are to: 

1 . Compile a detailed database of torques required to 
bend the joints of a space suit, using realistic, multi- 
joint human motions. 

2. Develop a mathematical model of the constitutive 
relations between space suit joint torques and joint 
angular positions, based on experimental data and 
compare other investigators’ physics-based models 
to experimental data. 

3. Estimate the work envelope of a space suited astro- 
naut, using the constitutive and compatibility rela- 
tions of the space suit. 

The body of work that makes up this report includes 
experimentation, empirical and physics-based modeling, 
and model applications. A detailed space suit joint 
torque-angle database was compiled with a novel exper- 
imental approach that used space-suited human test sub- 
jects to generate realistic, multi-joint motions and an 
instrumented robot to measure the torques required to 
accomplish these motions in a space suit. Based on the 
experimental data, a mathematical model is developed 
to predict joint torque from the joint angle history. Two 
physics-based models of pressurized fabric cylinder 
bending are compared to experimental data, yielding 
design insights. The mathematical model is applied to 
EVA operations in an inverse kinematic analysis cou- 
pled to the space suit model to calculate the volume in 
which space-suited astronauts can work with their 
hands, demonstrating that operational human factors 
metrics can be predicted from fundamental space suit 
information. 

2.0 Literature Review 

The extravehicular mobility unit (EMU) space suit, 
which is currently used for NASA EVA’s, has both hard 
fiberglass and soft fabric components. Mobility features, 
such as pleats that open as joints bend and rotational 
bearings, are built into all modem space suits. Without 
these mobility features, a person in a space suit would 
be virtually immobile. Even though space suits are 
designed to allow mobility, they restrict the wearer’s 

motion in significant and complicated ways. 20 


Limited space suit joint torque-angle data has been 
reported in the literature. Studies that used human sub- 
jects wearing space suits reported higher torques than 
studies that measured torques on joints of empty, pres- 
surized space suits. Higher torques with human subjects 
may be expected because contact between the space suit 
and the wearer’s body affects the deformed shape of the 
space suit, but it is unclear whether the discrepancy in 
experimental results is due to shortcomings in experi- 
mental methods or actual differences in observed 
torques. 2, 7 ’ 11 ' 28 

This research effort employed several mathematical and 
physical modeling techniques that were originally 
developed for other applications but are relevant to mod- 
eling an astronaut working in space, and specifically to 
provide an unique model of the torque-angle character- 
istics of space suit joints. Mathematical modeling tech- 
niques allow torques to be predicted from angle histories 
with accuracy, but, since they do not incorporate physi- 
cal principles, they do not contribute insights for design- 
ing space suit joints with better mobility characteristics. 
Three mathematical techniques for modeling hysteretic 
systems, the Krasnoselski-Pokrovski model 21 , the Prei- 
sach model, 26 and the Tao and Kokotovic model 35 , 
which were originally used for modeling magnetization 
and shape memory alloy actuators, were evaluated for 
possible use. The Preisach model was chosen because it 
can produce output curves similar in shape to the space 
suit torque vs. angle curves and its identification process 
is relatively simple. 

Physical models, in contrast to mathematical models, 
can lead to insights into the physical processes that gov- 
ern space suit mobility. Two physical models of the 
bending characteristics of inflated cylinders may be rele- 
vant to space suit joints. The beam model treats the 
inflated cylinder as a beam with a fabric wall that 
stretches, maintaining a constant internal volume. The 
mobility of the space suit joint is determined by elastic 
behavior of the fabric wall. 22 ' 25 The membrane model 
treats the fabric shell as an inextensible membrane. 
Bending deflections of the cylinder result in shape and 
volume changes. The work required to bend the joint is 
entirely due to compression of the gas inside the tube 
910 Differences in torque predictions between the beam 
model and the membrane model illustrate the relative 
importance of elasticity and gas compression in space 
suit joint mobility. 

One application of space suit joint mobility models is in 
calculating human factors performance metrics, includ- 
ing the one-handed space suited work envelope, which 
indicates the volume within which a space-suited astro- 
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naut can comfortably work with one hand. The current 
NASA space-suited work envelope 1 , which is used in 
planning EVA’s, differs qualitatively from other work 
envelopes found in the literature. 29 A recently-devel- 
oped computational technique allows suited work enve- 
lopes to be calculated using mathematical models of the 
torque-angle characteristics and kinematics of a space 


3.0 Space suit mobility database 

3.1 Experiment overview 

The objective of the experimental portion of this 
research effort was to obtain a quantitative database of 
joint angles and torques required to move a space suit’s 
joints, under realistic conditions. The impossibility of 
directly measuring joint torques in suited humans neces- 
sitated an indirect measurement approach, using both 
human subjects and an instrumented robot to obtain 
torque data. 

The human test subjects carried out arm and leg motions 
both wearing the space suit and not wearing the space 
suit, supplying realistic joint angle trajectories for each 
of 20 motions. The joint angle trajectories produced by 
the human subjects were then used as command inputs 
for the robot, so that the robot imitated the humans’ 
motions while torques were measured at each of the 
robot’s joints. Torques on the robot’s joints due to the 
weight of the robot and space suit, were subtracted, 
resulting in a consistent set of joint angle and torque 
data. 

3.2 Space suit description 

The space suit used in these experiments was an Extra- 
vehicular Mobility Unit (EMU) supplied under contract 
by Hamilton Sundstrand. A size large hard upper torso 
(HUT) was used.The space suit that was tested was as 
close as feasibly possible to a space suit that would be 
used in flight, however there were some differences 
between this space suit and a flight-qualified space suit. 
The space suit tested was designated as Class III hard- 
ware, approved for demonstrations or non-hazardous 
testing. Class III space suits are known to be less stiff 
than Class I, or flight qualified suits, because they are 
generally older and have been used more. Additionally, 
the scye bearings on the space suit used in these experi- 
ments, which connect the arms to the HUT, were 
mounted on a single-fold bellows, which allows some 
additional shoulder range of motion. This bellows con- 


figuration, known as a pivoted HUT, is no longer used in 
flight. 

Life support for the suited test subjects was provided by 
supplying breathing air on a non-recirculating vent loop. 
Air was supplied to the suit inlet by scuba tanks and 
exhausted from the suit to maintain a minimum air flow 
rate. One set of scuba tanks lasted approximately one 
hour. A volumetric mock-up was substituted for the por- 
table life support system (PLSS), resulting in an overall 
weight for the space suit and PLSS mock-up of 102 lb. 
(46.4 kg). In all cases, the space suit was pressurized to 

4.3 psi (30 kPa) above ambient pressure. 

3.3 Data collection with human subjects 

Subjects Four male subjects participated in the 
experiment; two of the subjects were experienced in 1 -g 
testing of the EMU. Subjects performed arm and leg 
motions in two sessions. The first session was conducted 
without the space suit and the second session was per- 
formed with subjects wearing the space suit, one week 
later. One of the inexperienced test subjects did not con- 
duct the simple motions described in Section for space- 
suited data collection because of a scheduling con- 
straint. The experiment was approved by MIT’s Com- 
mittee on the Use of Humans as Experimental Subjects 
and informed consent was obtained from each subject 
prior to each experimental session. 

Motions Human test subjects carried out ann and leg 
motions while wearing the space suit. Twenty motions 
w'ere used for human and robot data collection. The 
motions, listed in Table 1 , included 1 1 simple motions 
that isolated each of the robot’s degrees of freedom 
(wrist rotation was not used) and 9 complex, free-form 
motions, which included reaching over the head, across 
the body and to an object 50 cm above the floor. Sub- 
jects also walked on a treadmill and traversed a 12 cm 
high step to simulate planetary EVA tasks. Trial order- 
ing was randomized, with a first phase of simple 
motions followed by a second phase of complex 
motions. Subjects were instructed to move through the 
maximum range that they could comfortably reach. 

They were not instructed to attempt to achieve the maxi- 
mum range of motion. Free-form motions were included 
for two reasons: to obtain data on multi-joint, non-pla- 
nar motions and to compare the subjects’ choices of arm 
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and leg positions between suited and unsuited condi 
tions. 


for suited and unsuited configurations are shown in 
Figure .2 


TABLE 1. Arm and leg motions used in 
experiment 


Simple Motions 


Complex Motions 


Shoulder flexion 


Arm swing forward- 
backward 


Shoulder abduction 


Arm swing side to side 


Humerus rotation 


Leg swing forward back- 
ward 


Elbow flexion 


Leg swing side to side 


Hip flexion 


Overhead reach 


Hip abduction 


Thigh rotation 


Knee flexion 


Locomotion over 12 cm 
step 


Ankle rotation 


Locomotion on treadmill 


Ankle flexion 


FIGURE 2. Marker configurations for A.) suited and 
B.) unsuited configuration. 

The Cartesian positions of the markers were converted 
into joint angles for 1 1 arm and leg joints, which com- 
pletely describe the kinematics of the right arm, leg, and 
foot. Joint angles are measured with respect to the body 
segment closer to the torso or with respect to the torso in 
the case of the shoulder and hip joints. When all joint 
angles are zero, the arm and leg are both parallel to the 
torso and the palm of the hand and the toe both point 
forward. Each of the joint angles, with its zero position 
and positive direction, are shown in Figure 3. Figure A 
shows the five flexion angles: shoulder flexion (sf), 
elbow flexion (ef), hip flexion (hf), knee flexion (kf), 
and ankle flexion (af). Figure B shows shoulder abduc- 
tion (sa), hip abduction (ha), and ankle inversion (ai). 
Figure C illustrates the thigh rotation angle (tr), when 
hf=90 degrees and kf=90 degrees. Figure D and 
Figure E are both overhead views, showing, in Figure D, 


Ankle inversion 


Figure 1 shows an example of a subject performing the 
low reach trial. In the low reach trial, the subject’s task 
is to reach down and to the left, touching the corner of 
the step with his right hand. 


FIGURE 1. Low reach trial 


Data acquisition Kinematic data was collected on 
the subjects’ right arm and leg, using a Multitrax optical 
motion capture system (Adaptive Optics Associates, 
Cambridge, MA). Reflective markers were placed on the 
subject’s right arm, right leg and the mock-up PLSS, as 
shown in Figure . Outputs from the motion capture sys- 
tem were x, y, and z positions of each of the markers in a 
laboratory-based coordinate system. Marker positions 
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the humerus rotation angle (hr), with ef=90 degrees, and 
in Figure E, the ankle rotation angle. 



FIGURE 3. Joint angle conventions used for 

kinematic data for human test subjects 
and robot. 

3.4 Robot data collection 

Robotic Space Suit Tester, or M. Tallchief 

Known in our lab as M. Tallchief, a robotic space suit 
tester was built by Sarcos, Inc. (Salt Lake City, UT) for 
space suit mobility testing under the Small Business 
Innovative Research (SBIR) program for NASA 
Johnson Space Center. The robot, or M. Tallchief, was 
loaned to MIT by NASA in 1998 for space suit mobility 
research. 

M. Tallchief has 12 hydraulically actuated joints on the 
right arm and leg and 12 posable joints on the left arm 
and leg. At each actuated joint, potentiometers measure 
joint deflection and strain gauge load cells measure 
torque. Locations of the robot’s joints are indicated in 
Figure 4 and Table 2. 



FIGURE 4. The 12 degrees-of- freedom of M. 
Tallchief ’s joints. 


TABLE 2. RSST joint degrees of 
freedom 


Number 

Joint 

1 

shoulder flexion 

2 

shoulder abduction 

3 

humerus rotation 

4 

elbow flexion 

5 

wrist rotation 

6 

hip flexion 

7 

hip abduction 

8 

thigh rotation 

9 

knee flexion 

10 

ankle rotation 

11 

ankle flexion 

12 

ankle inversion 


The robot is suspended from a crane and supported at 
two points: a bolt at the “head” and a cable attached to 
the back, near the center of the torso. The torso is con- 
sidered to be the ground segment and remains in a fixed 
position, while the arm and leg move with respect to the 
torso. Adjustment of the cable supporting the torso 
allows the robot’s pitch angle to be adjusted and two 
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positioning joints within the torso allow forward/back- 
ward and left/right bending of the torso segment. 

The actuated joints of the robot are powered by an MTS 
model 510 hydraulic pump. Hydraulic fluid circulates in 
a loop from the pump, through each robot joint actuator, 
then back to the pump. All of the hydraulic lines and 
electrical cables exit the robot through a seal at the top 
so that a space suit may be pressurized on the robot. 

Control of robot motions is accomplished through ana- 
log circuitry and two PC computers, as shown in 
Figure 5. An EPC-5 486 computer supplies the user 
interface, and an EPC-6 386 computer performs the 
lower-level control functions. Analog control loops are 
closed on position, velocity, and torque on circuit boards 
in the Advanced Joint Controller (AJC) cage, with con- 
trol gains supplied by the EPC-6 through a digital inter- 
face. The AJC cage uses a modular architecture, with 1 2 
boards, one for each joint. The EPC-5 computer runs a 
Microsoft Windows 3.1 based user interface, called the 
Robotic Space Suit Tester Application (RSSTA), which 
passes commands to the robot at an update rate of 5 Hz 
and records position and torque data to disk. At the time 
the robot was built, the analog/digital command and 
control architecture allowed for faster loop closure than 
would have been possible by closing control loops in the 
digital domain in a computer. Although the analog posi- 
tion, velocity, and torque servo loops are hard-wired on 
the circuit boards, the computer interface does allow for 


time- varying control gains, a capability which has not 
been implemented in this study. 

The combination analog-digital control system is very 
robust to computer crashes, which happen periodically 
in Windows 3. 1 . If the EPC-6 computer stops receiving 
inputs from the EPC-5 computer, it holds the current 
commands to the robot’s joints. All of the control loops 
remain closed, so stability and performance of the 
robot's control system are preserved when the EPC-5 
computer crashes. In addition, the modular construction 
of the analog circuitry in the AJC cage aids in trouble- 
shooting. Boards for different joints may be exchanged 
for diagnostic purposes by changing jumper settings and 
inserting the boards in different slots in the AJC cage. 

The user interface program allows the robot to be com- 
manded in two ways: joints can be manually positioned 
by clicking on arrows in the positioning window or 
multi-joint trajectories can be loaded from files and exe- 
cuted. Joint positions and torques, sampled at 5 Hz, can 
be plotted on the screen and saved to files. Trajectories 
for the robot to follow may be created interactively, by 
manually positioning the robot and saving a series of 
trajectory points. The list of robot positions can be saved 
as a file and run at user-set speeds. Trajectories can also 
be created outside of the RSSTA application. 


Position 

Position Command Position 



FIGURE 5. RSST command and control. 

Robot operations 
Trajectory generation 

M. Tallchief was designed to approximate the joint axes 
and range of motion of a human arm and leg, but due to 
geometrical constraints in its construction, the robot has 
less range of motion in some joints than humans do. In 
particular, limits on two of the shoulder degrees of free- 
dom and one of the hip joints are quite different from 
those of human joints. The shoulder flexion joint does 
not allow negative shoulder flexion angles, which bring 


the arm past vertical, behind the body. The shoulder 
abduction joint also does not allow negative shoulder 
abduction angles, which bring the arm towards the body. 
Like the shoulder flexion joint, the hip flexion joint does 
not allow the hip to be flexed beyond vertical, behind the 
body. A complete list of the robot’s range of motion for 
each joint is given in Table 3, using the joint angle defi- 
nitions from Figure .6 

The robot’s limited range of motion prevents the raw 
kinematic data files collected from the human test sub- 
jects from being used directly to command the robot. 
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Commanding the robot with a joint angle that is beyond 
its range causes the joint to move to its limit and 
abruptly stop. To avoid spurious torque data caused by 
sudden changes in robot joint velocity, the data files col- 
lected from the human subjects were pre-processed so 


that all of the commanded joint angles were within the 
robot’s range of motion. 

TABLE 3. Robot joint angle limits 


Joint 

Maximum angle 
(deg) 

Minimum angle 
(deg) 

Shoulder flexion 

180 

-15 

Shoulder abduction 

90 

0 

Humerus rotation 

90 

-90 

Elbow flexion 

130 

0 

Hip abduction 

45 

0 

Hip flexion 

100 

0 

Thigh rotation 

22 

-22 

Knee flexion 

130 

0 

Ankle rotation 

20 

-20 

Ankle flexion 

30 

-45 

Ankle inversion/eversion 

20 

-20 


There are several methods that may be used to force 
commands to the robot to stay within the robot’s limits. 
The simplest method is to uniformly reduce the ampli- 
tude of the entire joint trajectory so that the maximum 
displacement falls within the limits. The problem with 
this approach is that it modifies the joint angle values for 
the entire data file, even though the vast majority ofjoint 
angle values are permissible. A more complicated 
scheme that minimizes changes to valid joint angles and 
discontinuous velocities while forcing the output joint 
angle trajectory to remain within the robot’s limits was 
implemented instead. 


f(t), and outputs y(t), the modified command, which 
does not exceed the robot’s limits for that joint. Two 
robot joint limit values were defined, as shown in 
Figure : b is the maximum possible robot joint displace- 
ment, and a is 90% of the way from the midpoint of the 
robot’s joint range to b. Values of the input, x(t), which 
exceed b are reduced to fall between a and b, while 
input values which are less than a are not changed. The 
maximum value of the input, x( t), is called x m . A scaling 
factor f(t) that accomplishes this transformation is given 
in Equation 1 , where k and c are constants whose values 
are set by matching boundary conditions. 


A nonlinear scaling method was used to eliminate robot 
joint commands that are outside the limits, while leaving 
non-extreme trajectory angles unchanged. This scheme 
forces the robot's commands to be in bounds by com- 
pressing the out-of-bounds angle values into a region 
near, but inside of, the joint angle limits. Compressing 
the out-of-bounds data rather than deleting it reduces 
velocity discontinuities, which may cause erroneous 
torque measurements. 


x(t)<a f(t) = 1 
x(t)>a f(t) = 


1 +k 


(x-c) 


Eq 1 


The constants k and c are set so that at x=a, y=x and at 
x=x m , y = b . The f(t) that satisfies these conditions is 
given by Equation 2. 


The pre-processing scheme takes a desired robot com- 
mand, x(t), as an input, multiplies it by a scaling factor 
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trates that y is always less than b and y is equal to the 
input value, .v, when x is less than a. 


yO) 

x(t) < a 
x(t)>a f(t) 


AOx(t) 

m = i 


. (*m- b ){x(t)-a) 

(x m -a) b 


Eq 2 


The results of applying this transformation are shown in 
Figure 6. For 0<t<tj, the input value, x, is less than a. 
The corresponding /f/K /■<?/,) is equal to 1 and conse- 
quently, the output y(0<t<t/)=x(0<t<f/).When x is 
greater than a, for tj<t<t 2 ,f is less than 1, and reaches a 
minimum of f=b/x m when x is at its maximum of x m . 
Thus, the maximum output value is equal to b, the maxi- 
mum allowable command. The plot oiy(t) vs. t illus- 


Equation 2 is valid only for non-zero, positive a and b. 
For negative or zero a and b, the same transformation 
can be accomplished by mirroring x(t) about the joint 
range midpoint, m, using Equation 3 to obtain x ’(t). The 
mirrored input x ’(t) is then transformed according to 
Equation 2, and the output, y ’(!) is then mirrored back 
according to Equation 3. 

x' = 2m -x y = 2m-/ Eq 3 



FIGURE 6. Trajectory motion limiting algorithm. 

The pre-processed commands for the robot are then for- 
matted into trajectory files that the RSSTA application 
can read. The robot’s trajectory files are ascii text files, 
which are normally created in the RSSTA application 
when the user manually positions the robot and saves 
trajectory points. Trajectory files were created from pre- 
processed data with the Matlab script trajconvert.m, 
which emulated the file format of the trajectory files cre- 
ated in RSSTA. Each trajectory file begins with three 
lines of text. Each line ends in a line feed character (ctrl- 
j or ascii 10). Following the three text lines are the com- 
manded joint positions. Each line contains joint posi- 
tions for the 12 joints, with up to 3 digits before the 
decimal point and 1 digit after it. Numerical values are 
aligned horizontally on decimal points. Numbers are 
separated by spaces so that there are always 6 characters 
between each successive decimal point, 4 characters 
before the first decimal point in the line, and 2 charac- 


ters after the last decimal point. Every line ends in a line 
feed (ctrl-j or ascii 10) and the end of file character is a 
carriage return (ctrl-m or ascii 13). If a trajectory file 
deviates from this format, either the RSSTA application 
will hang or it will load a short trajectory with incorrect 
joint angle values. 

Trajectory speed setting 

The speed at which the robot moves through a trajectory 
is set by the operator in RSSTA at the time the trajectory 
is executed, which allows the same trajectory to be run 
at various speeds. Speed setting is accomplished by 
entering a single velocity value in degrees per second, 
between 1 deg/sec and 50 deg/sec. The operator-entered 
velocity is assigned to the maximum velocity of the 
robot joint that has the largest angle range in that trajec- 
tory. For this experiment, the robot should move through 
the trajectory at the same speed as the human subject 
did. In order to accomplish this, the joint with the maxi- 
mum range and its maximum velocity would have to be 
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determined beforehand, stored, and correctly entered by 
the robot operator each time the trajectory is run. 

In order to reduce the likelihood of errors, a different 
method was used to set trajectory speeds. The wrist rota- 
tion joint mechanism had been removed from the robot, 
to allow the robot’s arm to fit in the EMU sleeve, as 
described in Section . A dummy command to the wrist 
rotation joint is included in all trajectory files for motion 
from -90 degrees to +90 degrees at a constant velocity of 
50 deg/sec. The wrist rotation joint then has the maxi- 
mum range and its velocity is set by the operator’s speed 
input. To ran trajectories at the same speed as the human 
subjects moved, the operator enters 50 deg/sec for all 
trajectories. 


Torque limits 

Torque limits for each joint are set and enforced by the 
RSSTA software. When a motion trajectory is running 
and one or more of the joint torques exceed the limits, 
the robot’s motion stops and the robot is commanded to 
return to the starting position of that trajectory. A dialog 
box is displayed indicating the joint that triggered the 
torque limit violation. Only one dialog box appears, 
even if multiple torque violations occur. Torque limits 
are not checked in manual operation. The joint torque 
limits used in the experiment are listed in Table 4. The 
limits used for wrist rotation are necessary because the 
wrist rotation load cell is disconnected and the open cir- 
cuit returns a constant value of -819 in*lb for wrist rota- 
tion torque. 


TABLE 4. Robot torque limits 


Joint 

Minimum Torque 
(in*lb) 

Minimum 
Torque (Nm) 

Maximum Torque 
(in*lb) 

Maximum 
Torque (Nm) 

Shoulder flexion 

-400 

-45.2 

350 

39.6 

Shoulder 

abduction 

-400 

• -45.2 

350 

39.6 

Humerus 

rotation 

-250 

-28.3 

250 

28.3 

Elbow flexion 

-350 

-39.6 

275 

31.0 

Wrist rotation 

-900 

-102 

500 

56.5 

Hip abduction 

-400 

-45.2 

2000 

226 

Hip flexion 

-1300 

-147 

1300 

147 

Thigh rotation 

-500 

-56.5 

500 

56.5 

Knee flexion 

-750 

-84.8 

750 

84.8 

Ankle rotation 

-500 

-56.5 

500 

56.5 

Ankle flexion 

-700 

-79.1 

700 

79.1 

Ankle inversion/ 
eversion 

-500 

-56.5 

500 

56.5 


Space suit installation space suit from damage. These included replacement of 

Although the robot had been designed for space suit structural parts in the robot and removal of non-essential 

mobility research, a space suit had never previously assemblies for ease in installing the space suit, 

been installed. Several modifications were made to the 
robot to allow the EMU to fit properly and protect the 
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The structural component between the top of the robot’s 
torso and the electrical-hydraulic feed through interface, 
referred to as the neck, was replaced to allow a gas-tight 
interface with the EMU space suit. The new neck has an 
angle of 25 degrees between its top and bottom surfaces, 
replacing the original neck which has a 55 degree angle 
between its top and bottom, as shown in Figure 7. 

The neck piece provides the mechanical interface 
between the EMU and the robot. When the space suit is 
installed on the robot, an aluminum plate, called the 
neck plug, is mounted to the neck ring of the HUT, 
where the helmet attaches. The aluminum plate is bolted 
to the robot’s neck piece, so that the space suit is sup- 
ported by the neck ring, rather than the inner surface of 
the HUT at the shoulders. The neck plug has a hole in 
the center for hydraulic hoses and electrical cables to 
pass through. Sealing is accomplished with an airtight 
seal between the EMU neck ring and the neck plug and 
a rubber gasket between the neck plug and the robot’s 
feed through interface, which is bolted directly above 
the neck plug. 


torso structure piece with a piece that was 3.2 cm 
shorter. The location of the modified torso piece is 
shown in Figure 8. A screw at the bottom of the robot’s 
torso was also replaced with a low-head version, to fur- 
ther shorten the torso. 



A. Original neck 



FIGURE 7. A.) Original neck with angle of 55 

degrees. B.) New neck with angle of 25 
degrees. 

In its original configuration, the robot’s torso was 
approximately 2.5 cm too long to allow the HUT and 
lower torso assembly (LTA) of the space suit to connect. 
The torso of the robot was shortened by replacing a 
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FIGURE 8. Modified torso structure piece, with 
original, larger torso structure piece at 
left. 

The robot’s lower arm also did not fit inside the space 
suit. The wrist rotation shaft and wrist disk were 
removed, shortening the lower arm length by 8.9 cm and 
reducing the arm end diameter from 7.6 to 6.4 cm. As a 
result, the wrist rotation joint was not functional, but it 
was never intended to be used in the experiment. 

In order to ease installation of the space suit on the 
robot, and eliminate the need to accurately size both 
arms and legs of the space suit, the robot’s non-actuated 
left arm was removed at the shoulder flexion pivot. The 
left foot was also removed at the ankle inversion pivot. 


the upper portion of the arm, which is normally horizon- 
tal, can be rotated to a vertical position so that the arm 
can be brought through the sleeve of the HUT. Care 
should be taken to ensure that the shoulder latch is 
engaged properly following installation of the HUT. The 
HUT was attached to the robot by bolts that pass 
through the robot’s feed through interface and the space 
suit’s neck plug into the robot’s neck piece. Using small 
lengths of threaded rod to align the bolt holes while the 
three parts are assembled aids this process. The robot 
was then reconnected to the crane, hydraulic lines, and 
electrical cables and raised to its normal position. The 
LTA, with the right boot removed, was pulled up over 
the legs and latched to the HUT at the body seal closure 
and the right boot was replaced. 


To protect the space suit from sharp edges and pinch 
points on the robot, a wet suit was placed on the robot 
before the space suit was installed. Foam padding was 
also used at the lower end of the plastic shell on the 
lower leg to prevent the space suit’s bladder material 
from being pinched between the foot shell and the lower 
leg during ankle flexion motion. A plastic cover was 
placed over the exposed end of the wrist rotation shaft. 


To install the space suit on the robot, the hydraulic lines 
and electrical cables were disconnected and the robot 
was disconnected from the crane. The robot’s feed 
through interface was then removed, down to the neck 
piece. With the robot held in a sitting position on the 
floor, the HUT was first installed on the robot. The arm 
was brought above the head and the shoulder latch was 
released by inserting an alien wrench through the access 
hole in the shell and pushing the latch mechanism, as 
shown in Figure 9. When the shoulder latch is released, 
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FIGURE 9. Shoulder latches. A.) Latch location. B.) 
Arm position with latch released. 

Position and torque data collection 

The EMU used in the human subjects phase of the 
experiment was installed on the robot and pressurized to 
30 kPa (4.3 psi), as shown in Figure 10. Air was sup- 
plied from a scuba tank to maintain the 30 kPa (4.3 psi) 
operating pressure in the space suit. Due to a high leak- 
age rate in the robot-space suit interface, one scuba tank 
lasted approximately 30 minutes. The high leakage rate 
required a suit depressurization and repressurization to 
replace the air tank every 30 minutes while data was 
being collected. 



FIGURE 10. EMU installed on robot. 


The space-suited robot performed 80 different motions 
that had been generated by the human test subjects. 
Robot joint angle and torque data for the 1 1 robot joints 
was obtained at both the speed that the human test sub- 
jects performed the motions and at half speed, to 
improve the robot’s trajectory-following performance. A 
total of 247 torque-angle data sets were obtained from 
the suited robot, including several replications of each 
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motion. In addition to the robot-generated angle and 
torque data, motion capture data was acquired for at 
least one replication of each motion, using the Adaptive 
Optics Associates optical motion capture system that 
had been used for data collection with human subjects. 

The robot also performed the above-listed motions with- 
out the space suit under two conditions: wearing the wet 
suit that had been used in space-suited data collection to 
protect the space suit from sharp edges on the robot and 
not wearing the wet suit. The unsuited torque and angle 
data was collected at full speed and half speed, the same 


FIGURE 11. Schematic of dafa collected from 
human subjects and robot. 

3.5 Data reduction 

The torque data that was recorded by the robot is a com- 
posite of torques required to support the weight of the 
robot’s limb, support the weight of the space suit arm or 
leg, and bend the joints of the space suit. Only the 
torques required to bend the space suit’s joints are rele- 
vant to the space suit model. Both experimental data and 
kinematic analysis were used to calculate the gravity- 


speeds at which the suited data had been collected. Over 
300 torque-angle data sets were collected from the robot 
without the space suit. Figure 1 1 shows in schematic 
form the nine motion capture and robot-generated data 
files that originate from one motion of the space-suited 
human subject. This report presents the space-suited 
human motion capture data and the robot-generated 
torque and angle data under suited and unsuited condi- 
tions. 


induced torques due to the robot and space suit. The 
gravity-induced torques were then subtracted from the 
joint torques recorded by the robot. 

Robot weight removal 

An empirical method was chosen for estimating torques 
induced by the robot's weight because the robot’s mass 
properties are not documented. Torques caused by the 
robot’s weight were measured by recording torque data 
while the robot performed all of the same motions that 
were used in suited robot data collection. The torques 
caused by the robot’s weight were eliminated from the 
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measured torque data by manually aligning the two files 
in time and subtracting the unsuited torque data from the 
suited torque data. Data collected with the robot wearing 
the wet suit only was used as the unsuited torque data. 
All data presented here has been corrected to remove 
torques induced by the weight of the space suit and 
robot. 

Space suit weight removal 

Torques caused by the space suit’s weight were esti- 
mated using kinematic analysis techniques, based on 
both the known mass properties of the space suit and the 
robot dimensions and joint axes. The kinematic analysis 
yielded the positions and orientation of each of the robot 
and space suit segments. Joint torques were calculated 
based the relative positions of the space suit segment 
centers of mass and the robot joints. 

The geometry of the torque calculation is shown in 
Figure 12. The vector from joint i to the center of mass 
of segment j is qj, q is the axis of joint i, nij is the mass 
of segment j, and g is the acceleration due to gravity. 
The torque caused by the weight of segment j about a 
point on z t is given by: 

I = tjj x mjg ■ Eq 4 

Only the torque along q is measured by the robot’s load 
cell at that joint, so the torque measured by the load cell 
on joint i is given by Equation 5. 

T U m bv xm jS)*?i Ec l 5 

The torque on each joint due to the suit weight is the 
sum of the 7jy torques for each segment that is further 
away from the torso than joint i is. 


each space suit segment, q, were calculated, giving the 
q and rjj vectors for a given set of joint angles. Coordi- 
nate axes used on each of the robot’s arm and leg seg- 
ments are shown in Figure 13. The Denavit-Hartenburg 
notation prescribes the locations of the x and z axes, 
based on the joint axes of the current segment and the 
previous segment. 5, 14 They axes are chosen to form a 
right-handed set of coordinates. 



FIGURE 12. Gravity-induced torque 7jy of segment j 
on joint i. 

Using the Denavit-Hartenburg notation for robot kine- 
matics, the position, q, and orientation, qof each of the 
robot’s joints and the position of the center of mass of 
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FIGURE 13. Coordinate axes for robot’s arm and 
leg. 

Based on the coordinate axes shown in Figure 13 and 
the joint angle definitions from Figure 3, the four 


Denavit-Hartenberg link parameters, a, d, a, and 6, for 
each robot segment are given in Table 5 for the leg and 
Table 6 for the arm. 


TABLE 5. Leg link parameters 


Link number 

fl(tn) 

d (m) 

a (deg) 

0(deg) 

1 

0.0699 

0 

-90 

hf 

2 

0 

0 

-90 

-ha-90 

3 

0 

0.40 

-90 

tr-90 

4 

0 

0 

90 

kf 

5 

0 

0.43 

90 

ar+90 

6 

0.0338 

0.0157 

90 

ai+90 

7 

0.0845 

0 

-90 

af 

8 

0 

0.220 

0 

0 
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TABLE 6. Arm link parameters 


Link number 

a (m) 

d (m) 

a (deg) 

0 (deg) 

1 

0 

0.127 

90 

sf 

2 

0 

0 

-90 

-sa+90 

3 

0 

0.318 

-90 

-hr-90 

4 

0.140 

0 

90 

ef-90 


The position and orientation of each segment are repre- 
sented by 4x4 homogeneous matrices. The homoge- 
neous matrix that transforms coordinate frame i-1 to 
coordinate frame i with rotation R and translation by 
vector x is given by: 

cosB. -cosu.sinB. sinu.sinB. a.cosB. 

i i i iii i 

sinBj. cosu f .costt. -sinaxosB^. a.sinki. 
0 sine i. cosa. d. 

° 0 0 1 

Positions and orientations are propagated down the 
kinematic chain by multiplying the ,4, /+t matrices. The 
homogeneous matrix representing the transformation 


from the base coordinate frame to the frame of segment 
n is given by Equation 6: 

n 

h, = Eq 6 

l 

By multiplying the A ,• ,+/ matrices corresponding to 
each segment, the position and orientation of each of the 
robot’s joints can be calculated. Positions of space suit 
segments are calculated by considering each space suit 
segment to be a shorter version of one of the robot’s seg- 
ments, with a length equal to the distance from the next 
robot joint closer to the torso to the space suit segment 
center of mass. To calculate a space suit center of mass 
position, the space suit segment is used as the terminal 
segment of the chain. Link parameters and masses of the 
space suit segments are given in Table 7. 


O' -i)< 


R x 
0 0 0 1 


TABLE 7. Space suit link parameters 


Space suit 
link 

Corresponding 
robot link 

£7(m) 

d( m) 

a(deg) 

^r(deg) 

mass (kg) 

upper arm 

3 

0 

0.225 

-90 

-hr-90 

0.677 

lower arm 
and glove 

4 

0.196 

0 

90 

ef-90 

2.25 

upper leg 

3 

0 

0.239 

-90 

tr-90 

0.741 

lower leg 

5 

0 

0.252 

90 

ar+90 

1.18 

boot 

8 

0 

0.0668 

0 

0 

2.24 


For each time step in the robot suited data file, the robot 
joint angle data is used to calculate the A, i+l matrix for 
each robot and space suit segment. The A,- j+ j matrices 
are multiplied to obtain the r, vectors that go from the 
origin of the base coordinate frame to the origin of the 
segment i coordinate frame, the z,- vectors that indicate 
the axis orientation of joint i, and the rj vectors that go 
from the origin of the base coordinate frame to the cen- 


ter of mass of space suit segment j. The robot and space 
suit link position and orientation vectors are substituted 
into Equation 5 to obtain the torque caused by the 
weight of space suit segment j on joint i. The space suit 
weight-induced torques are then summed for each joint 
to obtain the torque caused by the space suit's weight. 
Torques due to the space suit’s weight are then sub- 
tracted from the robot's torque data. 
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3.6 Experimental Results 
Example data 

After both the robot’s weight and the space suit’s 
weight are eliminated from the recorded torque data, the 
remaining torques indicate the torque necessary to bend 
the joints of the space suit. Examples of the space suit- 
induced torque plotted vs. joint angles are shown in 
Figure 14 - Figure 18. Each of the example plots repre- 
sents one test subject and one experimental trial, which 
includes several repetitions of the motion. Joint angles 
are measured according to the definitions provided in 
Figure 3. 


25 



FIGURE 14. Shoulder flexion torque vs. angle, 3 
repetitions. 



Shoulder abduction angle (deg) 


FIGURE 15. Shoulder abduction torque vs. angle, 3 
repetitions. 



Elbow flexion angle (deg) 


FIGURE 16. Elbow flexion torque vs. angle, 2.5 
repetitions. 

tflO- 



FIGURE 17. Hip abduction torque vs. angle, 3 
repetitions. 



Ankle flexion angle 


FIGURE 18. Ankle flexion torque vs. angle, 4 
repetitions. 

Database coverage 

The space suit database compiled in this work covers a 
larger range of joint angles than any other previous 
study on space suit joint torques. Angle ranges for which 
angle and torque data was obtained are shown in 
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Figure 19. The gray rectangles in Figure indicate the 
range of motion that was specified for the design of the 
EMU, or the robot’s range of motion for the humerus 
rotation, ankle rotation and ankle inversion joints, where 
the EMU range of motion was unspecified. The robot’s 
range of motion coincides with the EMU design specifi- 
cation in most joints, except the shoulder abduction 
joint. The robot’s shoulder abduction joint has a range of 
0 deg-90 deg, which is less than the EMU design speci- 
fication. The robot can, however, position its arm above 
its shoulder level by using the shoulder flexion joint, but 
it cannot measure true shoulder abduction joint torques 
in this configuration. Another instance in which the 
space suit design range of motion is much greater that 
the database range is the humerus rotation joint. The 


human test subjects did not exceed 50 deg of positive 
humerus rotation angle, even though the robot can 
accomplish a positive humerus rotation of 90 deg. Pub- 
lished data on the range of motion of unsuited individu- 
als for humerus rotation motions indicates that the 50th 
percentile humerus rotation limit in the positive direc- 
tion is approximately 60 deg, and the 95th percentile 

positive limit is 97 deg. 2 It is likely, then, that a 90 
degree humerus rotation would not have been achiev- 
able by the test subjects, even if they had not been wear- 
ing the space suit. 


sh flex 
sh abd 
elbflex 
hum rot 
hip flex 
hip abd 
knee flex 
thigh rot' 
ank flex 
ankrot 
ank Inv 





1 





-100 -50 ^ 0 50 100 150 200 

Suited angle range (deg) 


FIGURE 19. Angle range of space suit database 
compared to space suit design 
specification for joint range of motion. 

A total of 80 motion capture data sets were obtained 
from the space-suited test subjects, 247 torque data sets 
were obtained from the suited robot and 337 data sets 
obtained from the unsuited robot, resulting in over 580 
robot torque data sets. 


Error analysis 

Errors in the torque and angle values come from three 
sources: 

• Motion capture system errors in estimating the 
human subjects’s joint angles 

•Robot joint position errors due to imperfect track- 
ing of the trajectory 

•Robot torque measurement errors. 
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Motion capture errors 

Estimating the accuracy of off-the-shelf motion capture 
systems, such as the one used in this study, is difficult. 
Manufacturers make optimistic accuracy claims and suf- 
ficient computational detail is usually not available to 
the end user to make a detailed error analysis. This anal- 
ysis is an effort to make a reasonable, but somewhat 
conservative, estimate of the motion capture system’s 
accuracy in the configuration used in this study. 

The motion capture system measures the positions of 
reflective markers attached to the human test subjects’ 
arm, leg, and torso. Using the Cartesian positions of the 
markers, the angles between ann and leg segments are 
calculated. Random errors in measuring marker posi- 
tions can be propagated through the angle calculations 
to produce estimates of the angle errors that result from 
position measurement errors. 

The angle calculation is shown in Figure 20. The motion 
capture system outputs the x, y, z position of four points, 
A, B, C, and D, on two successive arm or leg segments. 
Using the positions of the four points, vectors rj and r 2 , 
along the two body segments, are calculated. The angle 
between these two body segment vectors, 8, is the joint 
angle. 



FIGURE 20. Joint angle calculation used for motion 
capture data. 

As shown in Figure 21, the vectors rj and r 2 are at 
angles 8j and d 2 with respect to an arbitrary coordinate 
reference. The joint angle 8 is the difference between 8j 
and & 2 . Focusing on a single body segment vector illus- 
trates how errors in measuring the positions of points A- 
D relate to angle errors. Figure shows body segment 


vector r h with point A’ displaced from point A by an 
error d, perpendicular to rj. 



A’ 


FIGURE 2 1 . Joint segment vector with measurement 
error 6 . 

The difference in angle between segment BA and seg- 
ment BA’ is given by 

60 * sin(50) = JiL = L*L Eq 7 

|d + § | Vu + 6 

If d 2 is the variance of the error in measuring r t , then 

*2 

var(60|) = — -. Eq 8 

r\ + b 

By the same reasoning, the variance of the error in 0 2 is 

&2 

var(bd 2 ) = -5 5 • Eq 9 

r\ + b 

Therefore, the variance of the error in measuring 8 is 
given by 

var(60) = var(60! ) + var(bQ 2 ) - -z 5 + T rEqlO 

r\ + b r 2 + b 

Equation 10 shows that errors in joint angle estimation 
increase when marker position estimation errors 
increase and decrease when the markers on a body seg- 
ment are separated by a greater distance. 

The manufacturer of the motion capture system claims 
that marker positions are determined with errors on the 
order of millimeters. The high accuracy claim may 
apply to the most favorable configurations, but the 
motion capture data acquired in this study is not likely 
to be that accurate. The markers used in this study were 
hemispheres of diameter 1 cm. It is reasonable to 
assume that, under good visibility conditions, the mark- 
ers are accurately located to within one diameter, or 1 
cm. 

Because the markers were blocked from the camera’s 
view at times, redundant markers were used. Each of the 
points A, B, C, and D in Figure is a position calculated 
by taking the means of the positions of two markers that 
were placed on a band around a test subject’s arm or leg. 
Averaging marker positions improves accuracy, but, 
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since markers were sometimes obscured, the full eight- 
marker accuracy was not obtained at all times. When 
two markers are available to be averaged, the variance of 
the error in calculating the position of point A, B, C or D 
is 0.5 cm , while when only one marker is available, the 
position error variance is 1 cm 2 . Table shows the angle 


error standard deviation for r t and r 2 dimensions typical 
of several joints, for different numbers of visible mark- 
ers. It is assumed that r } and r 2 span 0.75 of the length 
of the arm or leg segment. At least four markers must be 
visible in order to calculate an angle between the two 
segments. 


TABLE 8. Angle error estimates for different numbers of visible markers 



Joint angle error standard deviation (deg) 

Joint 

8 markers 

7 markers 

6 markers 

5 markers 

4 markers 

Elbow 

3.7 

4.2 

4.6 

5.0 

5.3 

Knee 

2.8 

3.3 

3.5 

3.8 

4.0 


Robot trajectory tracking errors 

The motion capture data from the human subjects was 
used as trajectories for the robot to follow in both suited 
and unsuited torque data collection. Errors in trajectory 
following were caused by loads imposed by the space 
suit and timing issues with the RSSTA software applica- 
tion. 

The space suit imposed significant loads on the robot as 
it attempted to follow the trajectories. As a result, the 
amplitudes of some of the motions were reduced on the 
robot, particularly in the shoulder abduction and hip 
flexion joints, where the loads due to the robot and space 
suit’s weight and the space suit’s stiffness were largest. 
The lack of stiffness of the robot’s joints effectively lim- 
ited the range over which the joint could be positioned 
when the space suit was on the robot. Because the 
RSSTA application does not allow the user to command 
joint deflections that are outside the robot’s limits, the 
actual joint angle that was obtained when the maximum 
joint angle was commanded was the effective limit. This 
problem can be remedied in the future by tuning the 
servo gains on the robot’s shoulder and hip joints to 
have higher position error gains. Effects of tracking 
errors when the robot was wearing the space suit are 
shown in Figure 22, where the percentage of unsuited 
joint angle range that was achieved in suited data collec- 
tion is plotted for each of the 1 1 simple motions. 



FIGURE 22. Percentage of unsuited range achieved 
in suited robot trials. 


An additional source of error in robot trajectory follow- 
ing arose from timing issues. When the robot was pro- 
grammed to follow motion capture data trajectories at 
full speed, the joint angle data recorded by the robot 
indicated in many cases that the robot actually moved 
slower than the command trajectory. The duration of the 
data file exceeded the duration of the trajectory file and 
it appeared that samples were being dropped while the 
robot was performing motions. The robot would hold all 
joints in the same positions for 2-3 samples at a time, 
then proceed with the motion. The number of “held” 
samples corresponded to the excess time taken in exe- 
cuting the motions. Less sample dropping occurred 
when trajectories were run at lower speeds. To produce 
better matching of the robot’s motions to the command 
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trajectories, suited and unsuited data was collected for 
all trajectories at both full speed and half speed. 

The angle data in the space suit torque-angle database is 
the robot’s actual position, not the commanded position, 
so tracking errors have the effect of slightly distorting 
the motions that the robot performed; tracking errors do 
not cause errors in the measured angles in the space suit 
database. Figure 23 shows root mean square (RMS) dif- 
ferences between suited and unsuited robot joint angles, 
plotted by trial type. The RMS errors, which range from 
0.5 deg to 3 deg, are of the same order of magnitude as 
the motion capture data errors listed in Table 8. 



Motion Ivoe 

FIGURE 23. Mean RMS angle differences between 
suited and unsuited robot motions. 

Torque measurement errors 

Errors in torque measurement arise from bias, noise, and 
quantization. The robot’s load cells were calibrated 
when it was originally set up at Johnson Space Center. 
These calibration factors were in use when suited and 
unsuited data was collected with the robot. The load cell 
calibrations were again evaluated after the data collec- 
tion by applying known torques to each load cell and 
comparing the RSSTA torque output to the known 
applied torque. The later assessment indicated that cali- 
bration factor errors were all less than 4%. 13 

Noise and quantization effects interact to result in a ran- 
dom error of approximately 0.1 13 Nm (1 inch pound) in 
all joints. Torque values are represented in the RSSTA 
software as integers in the units of inch pounds, reduc- 
ing the resolution by as much as a factor of 10, com- 


pared to the hardware-imposed resolution limits, in 
some cases. As a result, the resolution of the torque data 
depends on the calibration factor for each joint. The res- 
olution is the larger of the torque corresponding to 1 A/ 
D count or 1 inch pound (0.113 Nm). For all joints 
except the hip flexion joint, the torque resolution is 
0. 113 Nm; the torque resolution for the hip flexion joint 
is 0.226 Nm. 

3.7 Database contributions 

The objective of the experimental portion of this 
research effort was to compile an extensive database of 
the torques required to bend the joints of a space suit, in 
realistic, human-generated motions. The experiments 
were carried out in two phases, with both human test 
subjects and an instrumented robot. Human test subjects 
wearing the space suit performed 20 simple and com- 
plex motions, while joint angles of their arm and leg 
were recorded using video motion capture. The space 
suit was then installed on an instrumented robot, which 
was driven to reproduce the human subjects motions 
while joint torques were recorded. Because the robot 
had never been used with a space suit, several modifica- 
tions were made to the robot to allow the space suit to fit 
on the robot and protect both from damage. Contribu- 
tions due to the weight of the space suit and robot were 
subtracted from the torque data, resulting in a consistent 
set of space suit joint torques and joint angles in realistic 
motions, with angles accurate to approximately 2 deg-5 
deg and torques accurate to approximately 0. 1 Nm. 

The database compiled in this work is more extensive 
than any other published space suit torque-angle data- 
base, covering 1 1 joints over a large range of angles. 
Realistic, three-dimensional human-generated motions 
were used for data collection. Because torque data was 
collected by the robot as a surrogate for a human occu- 
pant of the space suit, the torques measured in this study 
are more representative of realistic conditions than data 
from previous studies that measured joint stiffnesses for 
empty, pressurized space suits. The space suit torque- 
angle database serves as a basis for developing and vali- 
dating both mathematical and physical models of space 
suit joint mobility characteristics. 

4.0 Modeling 

The space suit torque-angle database described in 
Section 3.0 presents a unique opportunity to develop 
models of space suit mobility and verify them against 
experimental data. Models of space suit mobility are 
useful in two applications: numerically predicting the 
torque required to bend a space suit’s joints and under- 


Quantifying Astronaut Tasks: Robotic Technology and Future Space Suit Design - PI: Dava Newman 


21 



standing the physical processes that determine how 
mobile a space suit’s joints are. These two applications 
require two different modeling approaches: first, a 
descriptive mathematical modeling technique based on 
experimental data and second, a theoretical model based 
on physical principles. Section 4.0 describes two model- 
ing efforts: a mathematical model based on empirical 
data that predicts the torque required to bend the space 
suit’s joints and a comparison between two physics- 
based models of bending pressurized cylinders and the 
space suit torque angle database compiled in 
Section 3.0. The opportunity to validate the mathemati- 
cal and physical models against experimental data is a 
unique aspect of this research effort. 

4.1 Mathematical modeling 

The Preisach model The Preisach hysteresis model 
reproduces a hysteresis curve by summing contributions 
from the simplest possible hysteresis transducers. The 
primitive hysteresis transducer, y(a, fi) , has an output 
value equal to either +1 or -1. For increasing inputs, the 
output switches from - 1 to +1 at an input value of a, and 
for decreasing inputs, the output switches from +1 to -1 
at an output value of ft. 

To construct more complicated hysteresis transducers 
with continuous, non-unity outputs, the Preisach model 
uses a weighted sum of simple hysteresis operators. The 
weighting function ju(a,/3) is defined as a function of the 
combination of upward and downward switching values, 
a, ft, of the hysteresis transducer. The Preisach func- 
tion, jtt(a,(3), is defined for all a > p, -a 0 < a < a 0 , and - 
a 0 < P < «o, forming a triangle in a- (5 space. 

Construction of the output of the composite hysteresis 
transducer is done by integrating the individual y„ p 
values, as shown in Equation 1 1 . 

/(0=JJ g(a, P)y K pu(/)tfcu/p Eqll 

(<x>P) 

Calculation of f(t ) is aided by a graphical representation 
of the a-p space. The weighting function fi(a,p) is 
defined over the triangle that is bounded by the a=fi line 
and the maximum values of a and ft, a 0 and p 0 , which 
are defined by the saturation limits of the output. To 
obtain f(t). Equation 1 1 is integrated over this triangle. 

To perform the integration in Equation 1 1 ,the triangle 
may be subdivided into two sets: 

•Region S+, where corresponding y u operators 
are “up”, or equal to +1 


• Region S-, where corresponding y (l p operators 
are “down”, or equal to -1 

Substituting +1 for y (t p in the S + region and -1 for 
y H p in the S’ region, Equation 1 1 simplifies to 

/(/) = J'JVfct, pjt/adp - J’J’gfa, $)dad$ . Eq 12 
s 

Equation 12 can be integrated if the boundary between 
the S~ and S’ regions is known. The input history is 
stored by drawing the boundary between the S + and S’ 
regions. The boundary is drawn by constructing line 
segments in a staircase pattern within the triangle based 
on the value of the input u(t) and whether the input is 
increasing or decreasing. The boundary is drawn 
according to the following rules: 

1 . The boundary starts on the a=a (J segment if the ini- 
tial input is descending and the P=-a 0 segment if 
the initial input is ascending. 

2. Subsequent boundary segments are drawn horizon- 
tally or vertically depending on whether u(t) is 
increasing or decreasing: 

• Increasing u(t): horizontal line segment at a=u 

• Decreasing u(t): vertical line segment at /3=u 

3. A boundary segment is obsolete if its a value is less 
than the a value of a later segment having the same ft 
value or if its p value is less than the ft value of a 
later segment that has the same a value. 

4. The last line segment ends on the a=p line. 

An exploration of the Preisach hysteresis model focused 
on four of the important features of the outer hysteresis 
curve shape: number of loops, direction of loops, sepa- 
ration between increasing and decreasing curves and dT/ 
du, resulting in several important insights. The proper 
choice of n(a,P) can produce single or multiple hystere- 
sis loops, in either clockwise or counterclockwise direc- 
tions. The separation between increasing and decreasing 
curves can be written as an integral of n(a,p) over a 
region in the top left comer of the a-fi space. Thus, high 
amplitudes of n(a,P) in the top left comer of the a-p 
space, where a is near a 0 and p is near -a 0 , result in 
large separations between the increasing and decreasing 
output curves. 

Model identification The model identification pro- 
cess for the Preisach model involves calculating the 
weighting function n(a,P) from experimental data. The 
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central premise of the Preisach model identification 
methods is that the difference between two output val- 
ues is equal to the integral of /j( a,/3 ) over a region whose 
bounds are known from the input history. Calculating 
the integrated for a sufficient number of regions 

allows the Preisach model to be implemented to predict 
the hysteretic system’s output as a function of its input 
and input history. The following derivation shows how 
the integral of n(a,fi) over a region with known bounds 
can be calculated from the difference between two out- 
put values. 

The input to the hysteresis transducer begins below the 
low input limit, -a 0 , and increases to a value u /, which 
is below the high input limit. The corresponding output, 
when u=uj, is 7=7). In the a-fi plane, the S + /S" bound- 
ary lies at a=U]. The input then decreases from u=Uj to 
u=u 2 , and a vertical segment is added to the S + /S" 
boundary at p=u 2 . The output when u=u 2 is 7=7). 

According to 12, the output T(u) is equal to the integral 
of n(a,P) over the S + region minus the integral of n(a,/3) 
over the S' region. Using 12 to calculate the difference 
between T ) and 7) results in the following expression 

for the output difference, in terms of the S + and S' 
regions at u=u/ and u-u 2 . 


explained by Ge and Jouaneh 11 , is based on the same 
premise, but avoids differentiating data in the identifica- 
tion step and double integrations in the implementation 
step. This method avoids differentiation and integration 
by calculating the integral of n(a,j3) over a collection of 
triangles in a-fi space from output differences, then 
uses sums and differences of the triangle integrals to 
construct the output for any input history. 

According to Equation 13, the integral of n(a,fi) over a 
triangle bounded by the a=fi line and the a=u /, /3=u 2 
point is equal to 7)- 7), as long as the input increased 
from its lower limit to uj, then reversed direction and 
decreased to u 2 , with no other direction reversals. This 
relationship allows the integral of n(a,fi) to be calcu- 
lated based on the appropriate output differences for tri- 
angles bounded by the a=fi line and any point inside the 
large triangle over which n(a,(i) is defined. The quantity 
X(a t ,Pi) is defined as the integral of n(ct,/3) over a trian- 
gle bounded by a,, P/ and the a=/S line. If w increases 
from its low limit, reverses direction at m=cq , then 
decreases to ^ { ,X(a h li,) is given by 

*(«,,{*,) = (T(u = a,)- 7’(u = p,)) 

“ lU| Eq 14 

*(a„P,) = 2jjg(a,p)<fa45 

Pi P 


7) -T 2 = T max -2 ffn(a, f})c/a</ff 
s' i 



max 


-2 J J|x(a, P )dad$ 

S' 2 


T\~ T i = 2 j / u(a, p)c/a<7p 

r 2 -s' , 


Eq 13 


The difference between integration areas, S^-S), for 
u=u 2 and u=uj, is a triangle that has vertices at (u i,u 2 ), 
(uj, Uj) and (u 2 , u 2 ). When the input decreases from uj 
to u 2 , the difference in output values is equal to the inte- 
gral of i-i(a,l3) over the triangle. 

Numerical Implementation To determine u(a,(5) 

from output differences, Mayergoyz 26 suggests differen- 
tiating the output data twice with respect to the input. 
Taking two derivatives of measured data values would 
amplify random noise to unacceptable levels, so the 
method that Mayergoyz 26 recommends is not practical. 

An alternative Preisach model identification scheme, 
developed by Doong and Mayergoyz 8 and further 


lfX(a,fi) is known for all -a 0 <a<a 0 , -a 0 <f5<a, 0} a>/3, 
then X(a,/3) values can be added and subtracted to con- 
struct any Preisach model output, provided that the 
boundary between S + and S' is drawn according to the 
rules described above. The integral of p(ct,P) over S + is 
given by: 

n 

T = -X(a 0 ,-a 0 )+ 2-l" + l X(a,,p,.) Eq 15 

/= i 

Error' Analysis Because the hysteresis model coeffi- 
cients are determined from experimental data, random 
errors in the experimental data lead to random errors in 
the model output whose statistical properties can be pre- 
dicted. The error analysis provides a method for gener- 
ating confidence intervals for the Preisach model output 
if the model is identified and implemented according to 
the Doong and Mayergoyz 8 method. 

According to Equation 1 5, the Preisach model output is 
the sum of positive and negative X(a,f3) values at the 
vertices of the S + /S' contour. Assuming that the errors in 
the experimental torque and angle data that was used to 
calculate X(a,(5) are random, white noise, errors in 


Quantifying Astronaut Tasks: Robotic Technology and Future Space Suit Design - PI: Dava Newman 


23 


X(a,f5) should be uncorrelated with errors in X at other 
values of a and (3. Thus, the variance of the model out- 
put T is equal to the sum of the variances of the individ- 
ual X( a, (5) values that were summed to obtain T. 

Summing the variances of the X values at the n S + /S' 
boundary vertices results in the variance of the model 
output T. The variance of the measured torque is err- 
and the variance of the measured angle is ct A . 

2 2 ^ (dX: dX \ 2 

var(T) = (2 n)a T + a A ^ + Eq 16 

4.2 Results: Hysteresis modeling 

Preisach model coefficients X( a,f5) were fitted to the 
first-order transition curves that were produced for space 
suit torque-angle data. Model coefficients were obtained 
for elbow flexion, hip abduction, hip flexion, knee flex- 
ion, ankle rotation, and ankle flexion motions. These 
experimentally-determined model coefficients were then 
used to predict the torques required to bend the space 
suit joints for motions produced by space-suited human 
subjects, then the torque predictions generated by the 
model were compared to experimental data. 

The Preisach model coefficients are defined in terms of 
the upward and downward switching thresholds, a and 
15, of the primitive hysteresis transducers. The model 
coefficients, X(a,f5), can be plotted in three dimensions 
vs. a and jfl, over the triangular region for which a and f5 
are defined: -a 0 <a<a 0 , -a 0 <(5<a 0 , a>/5. The experi- 
mentally-determined model coefficients for the elbow 
flexion, hip abduction, hip flexion, knee flexion, ankle 


rotation, and ankle flexion joints are plotted in Figure 24 
- Figure 29. 



Elbow flexion 


150 


150 0 


Beta (deg) 


Alpha (deg) 


FIGURE 24. Preisach model coefficients plotted vs. 
a and (3 for the elbow flexion joint. 

Hip abduction 



FIGURE 25. Preisach model coefficients plotted vs. 
a and p for the hip abduction joint. 
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Hip flexion 


Ankle rotation 



FIGURE 26. Preisach model coefficients plotted vs. 
a and (5 for the hip flexion joint. 



Knee flexion 


x 


FIGURE 27. Preisach model coefficients plotted vs. 
a and (i for the knee flexion joint. 
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FIGURE 28. Preisach model coefficients plotted vs. 
a and (I for the ankle rotation joint. 


Ankle flexion 



20 -20 


Beta (deg) 


Alpha (deg) 


FIGURE 29. Preisach model coefficients plotted vs. 
a and for the ankle flexion joint. 


The hysteresis models were implemented for each joint 
and compared to experimental torque-angle data for 
motions that were generated by space-suited human 
subjects. Joint angle data from the human subjects was 
used as an input for the hysteresis model, which gener- 
ated a torque prediction as well as a 95% confidence 
interval on the prediction. The experimental joint torque 
data is compared to the model predictions and confi- 
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dence intervals for the elbow in Figure 30 and for the 
knee in Figure 31. 



Time (sec) 

FIGURE 30. Elbow flexion torque compared to 

hysteresis model prediction for subjects 
B (top), C (middle), and E (bottom). 



o 5 10 15 

Time (sec) 


FIGURE 3 1 . Knee flexion torque compared to 

hysteresis model prediction for subjects 
B (top), C (middle), and E (bottom). 


The Preisach model outputs match well with data for 
elbow flexion, hip abduction and knee flexion, with 
R >0.6 for these joints. The error estimates generated 
appear to be a reliable, or at least conservative, estimate 
of the emor standard deviations. 


The Preisach hysteresis model is a useful tool for pre- 
dicting the torque needed to bend space suit joints as a 
function of time for complicated angle trajectories. The 
ability to reproduce realistic space suit torques dynami- 
cally makes the Preisach model well-suited for use in 
dynamic simulation of EVA tasks. The experimental 
data reported in Chapter 3 as well as other investigators’ 


data indicate that the torque-angle relationship for space 
suit joints is markedly hysteretic. Consequently, repro- 
ducing dynamic behavior of space suit joints requires a 
model that captures hysteresis; other modeling tech- 
niques such as linear regressions on angle and angular 
velocity will not accurately predict the effects of the 
angle history on the torque required to bend the space 
suit’s joints. 

The Preisach hysteresis model has some limitations in 
modeling space suit joints, however. To identify the 
model coefficients, it is necessary to vary the input 
between several distinct minima and maxima. Because 
of this, the model identification process is not well- 
suited to human-generated motions, although character- 
istics of human-generated motions, such as range and 
speed, can be incorporated into the input data for real- 
ism. Another limitation of the Preisach model is that the 
maximum input and output are set when the model coef- 
ficients are identified. When the model is implemented 
later, if the input exceeds the previously-set maximum, 
the output saturates and errors between the model pre- 
diction and actual torques may become large. Modeling 
the hysteretic characteristics of space suit joints is 
essential for making accurate predictions of the torques 
required to bend the space suit joints in dynamic situa- 
tions. This work demonstrates that the Preisach hystere- 
sis model accurately reproduces torques needed to bend 
space suit joints. 

4.3 Physics-based modeling 

Beam model The beam model, which was developed 
by Comer and Levy 3 and extended by Main, Peterson, 

and Strauss 22 ' 25 , treats a pressurized cylinder as a long, 
slender member, loaded in a single plane, whose behav- 
ior is governed by elasticity and buckling phenomena. 
The central idea of the beam model is that fabric can 
sustain only extensional stresses; when the fabric wall of 
the beam is compressed, it wrinkles and does not con- 
tribute to the stiffness of the beam. The beam model pre- 
dicts the extent of the wrinkled portion of the beam and 
integrates the extensional stress resultants over the ten- 
sioned, un wrinkled, portion of the beam to obtain a 
moment-curvature relationship for the pressurized 
beam. 3, 22 ‘ 25 Stress resultants are normalized by the 
thickness of the fabric, and expressed in terms of load 
per unit length in units of N/m, because fabric stress- 
strain behavior is not highly correlated with fabric thick- 
ness. 

The applied loads and the constitutive relations of fab- 
ric, which relate stresses and strains and are given in 
Equation 17, determine the extent of wrinkling in the 
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fabric. In general, fabric is assumed to be orthotropic, 
with Poisson ratio vand different moduli, E h and E l in 
the hoop and longitudinal directions. 


e 


L 




Eq 17 


The beam model predicts that the fabric wall of the 
beam wrinkles when the local longitudinal strain, e L , 
becomes negative. Equation 1 7 shows that the cross- 
coupling component in the constitutive relations makes 
the fabric wrinkle at positive values of longitudinal 
stress, a L , when the hoop stress, a H , is positive. The 
hoop loading due to pressurization thus makes the fabric 
wrinkle at lower applied loads as the beam is bent. The 
longitudinal stress at which wrinkling occurs is given 
by: 

a L = vpr Eq 18 

When the fabric is wrinkled, its stress resultant is con- 
sidered to be zero, because it cannot resist compressive 
loads. The longitudinal stress is maximum on the out- 
side of the bend and decreases linearly until it reaches 
the wrinkling limit, then the wrinkling condition sets it 
equal to zero. 

To determine the curvature of the beam, it is first neces- 
sary to find the extent of the wrinkled region, using the 
relation between applied moment, M, internal pressure, 
p, and beam radius. A polar coordinate frame is defined 
with 6=0 on the inside of the bend, at the center of the 
wrinkled region and 6=6 0 at the edge of the wrinkled 
region. The value of 6 0 is given by Equation 19. 

M ^[(it-6 0 )+ sine o cos0 o ] + 

sine 0 + (jt-e 0 )cose 0 Eq ^ 

-o[(jt - 0 O ) 2 - (re 0 o )sin0 Q cos0 o - (2sin0 o ) 2 ] 
sin0 Q + (ji - 0 Q )cos0 o 


The 6q that is calculated by solving Equation 19 numer- 
ically for specified M/pr 3 can then be substituted into 
Equation 20 to obtain the curvature resulting from the 
applied moment. 


A/-2i7?r 3 sin0 o 
Er [(jt-0 o )+ sin0 o cos0 o ] 


Eq 20 


In contrast to the beam model, the membrane model 
treats the fabric cylinder wall as an inextensible material 
that transmits forces only along its surface and only in 
tension. Bending deflections of the structure result in 
changes in its cylinder’s cross-sectional shape and the 


volume that it encloses. The cross-sectional shape of the 
fabric tube is detennined by the assumptions of inexten- 
sibility and exclusively tensile loading. These assump- 
tions result in reliable approximations of the tube’s 
shape, even when deflections are large. 

Membrane model The membrane model uses a 
variational principle to relate applied bending moment, 
force applied at the tube’s ends, and the bending angle 
of the pressurized fabric tube. The potential energy of a 
tube with volume V, internal pressure p, end force Q, 
bending moment M, bending angle <p and linear dis- 
placement <5 is given by: 

n = -pV-Qb-M^ Eq 21 

When the system is at equilibrium, the potential energy 
is minimized. For the case of pure bending, Q=0 and 
8=0. The moment-angle relationship that minimizes the 
potential energy is calculated by holding the bending 
moment, M, fixed and differentiating the potential 
energy, 77, with respect to the bending angle <p. When 

= 0 , the potential energy is minimized. 


dU 
d( j> 



-M = 0 


=> M 



Eq 22 


Equation 22 shows that the equilibrium bending 
moment-bending angle relationship is set by the rela- 
tionship between the internal volume of the tube and the 
bending angle. The assumptions of inextensibility and 
tensile loading set the shape of the tube as a function of 
bending angle. Based on the four principles of mem- 
brane shape listed in Reference 10 , the shape of the bent 
tube can be determined. The geometry of a bent pressur- 
ized tube is shown in Figure . For a tube of radius R, the 
radius of the bent outer surface of the tube is 2R. The 
inner surface at the bend folds inward, towards the cen- 
ter of the tube, forming a kink in the tube wall. Away 
from the kinked region, the tube is straight. Points P and 
S, which lie on the tube wall, are defined to aid in calcu- 
lating the tube’s cross-sectional area. Point P is located 
on the kink at the inside of the bend, and point S is the 
corresponding point on the outside of the bend. Point P 
and point S are the same distance away from the end of 
the tube. 6 is the angle between a line drawn from the 
inside of the bend to point S and another line that is nor- 
mal to the straight sides of the tube. The length of the 
bent portion of the tube on the outer edge is 4R<p. 
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FIGURE 32. Bent tube shape, according to 
membrane model. 

The locations of points P and S are defined as follows. 
x p = 0 

v P = 2R( 1 + 0) 

Eq 23 

x s = 2R(\ - cos(<j> — 0 )) 
y s = 2/f sin(cj) - 0) 

The distance between points P and S, which can be cal- 
culated from Equation 23, and continuity of slope 
between the rectangular portion and the circular portion 
of the cross-section allow the cross-sectional area to be 
calculated as a function of 6. The cross-sectional area of 
the tube is given by Equation 24, as a function of H, the 
distance between points P and S. 

A(Q) = Eq 24 

The cross-sectional area is integrated over the length of 
the tube to obtain the bent tube volume. The moment- 
bending angle relation is obtained by differentiating the 
volume numerically with respect to <f>, according to 22. 

A comparison of the torque-angle predictions of the 
beam model, membrane model and experimental data 
from the knee joint is shown in Figure . The membrane 
model agrees with experimental data within 30 degrees 
of the knee joint’s equilibrium angle, while the beam 


model’s torque predictions do not agree with the experi- 
mental data. 



FIGURE 33. Comparison between beam model, 

membrane model and experimental data 
for knee joint. 


The beam model and membrane model represent oppo- 
site extremes of the processes governing space suit joint 
mobility. The beam model assumes that the space suit 
joint torque-angle relationship is determined solely by 
elastic defonnations of the fabric wall of the space suit, 
with no changes in the internal volume of the space suit 
segment. The membrane model makes the opposite 
assumptions: that the space suit fabric does not stretch, 
although the cross-sectional shape deforms, resulting in 
reduction in internal volume as the joint is bent. Accord- 
ing to the beam model, bending moments only stretch 
the space suit fabric; according to the membrane model, 
bending moments only compress the gas inside the 
space suit. 

It is reasonable to expect that the actual behavior of a 
space suit falls between these two extremes, that both 
elasticity and volume changes determine space suit joint 
mobility. Comparing these two approximate models to 
experimental data illustrates where on the continuum 
between elasticity and volume change the space suit 
joint behavior falls and indicates whether the space suit 
behavior can be approximated by one model or the 
other. 


The beam model is inconsistent with experimental data 
for the elbow and knee joints, both in the bending 
moment magnitudes predicted and the trend in bending 
moment with increasing deflection. The beam model 
predicts bending moments for the elbow and knee joints 
that exceed the experimental data by at least a factor of 
5. In addition, the beam model predicts a decreasing 
slope of the moment-angle curve as bending angle 
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increases, while the data shows a increasing slope with 
increasing bending angle. 

The shape of the torque-angle curves generated by the 
membrane model is consistent with the experimental 
data, because the model’s torque-angle curves are flat 
near the equilibrium position of the joint and increase in 
slope as deflection increases. The membrane model fol- 
lows the experimental data over a bending angle range 
approximately 30-50 degrees from the equilibrium 
angle, for the elbow and knee joints. Outside of this 
range, the amplitude of the experimental torque data is 
greater than the membrane model’s predictions. 
Although the membrane model can predict the torque 
needed to bend the joints over about half of the available 
range of motion, the model’s underprediction of torques 
near the extremes of the angle range prevent it from 
being used to predict a joint’s angular range of motion. 
Other physical processes, including fabric bunching, 
sideways compression of the space suit ann or leg, and 
friction between fabric layers, which were not included 
in the membrane model, likely play an important role in 
determining the torque needed to bend the space suit 
joints at high bending angles. 

The results of this analysis indicate that elastic deforma- 
tions of the space suit fabric do not contribute signifi- 
cantly to the torque needed to bend the elbow and knee 
joints near the equilibrium angle, and the mobility of the 
EMU elbow and knee joints near their equilibrium 
angles is not sensitive to changes in the material proper- 
ties of the fabric. The membrane model predicts that 
only internal pressurization and joint geometry deter- 
mine mobility. 

5.0 The Work Envelope: Applying 
space suit modeling to EVA 
operations 

An important issue in planning extravehicular activity 
(EVA) tasks and evaluating EVA worksites is determin- 
ing whether astronauts can reach and comfortably work 
at the desired worksite. A computational work envelope 
analysis was developed, which utilizes experimental 
joint angle and torque data measured during space suit 
experiments and inverse kinematics analysis, to deter- 
mine the volume in which a space-suited astronaut can 
comfortably work with one hand. For a point to be 
within the workspace, it must be visible and not require 
excessive torques at the shoulder or elbow joints. These 
constraints were applied by calculating the joint angles 
required to place the hand on a designated target, deter- 
mining the torques at the shoulder and elbow that were 


needed to hold the hand position, then applying joint 
torque and visibility limits. Work envelopes were gener- 
ated from both male and female anthropometric data. A 
sensitivity analysis on joint range of motion and visibil- 
ity limits showed that increasing the range of the shoul- 
der abduction/adduction joint to allow the arm to move 
closer to the body centerline greatly increases the work 
envelope volume. Increasing visibility in the upward 
and downward directions leads to sizable increases in 
work envelope volume, while visibility increases in the 
left and right directions do not increase the work enve- 
lope size. 

A critical issue that is addressed in planning for extrave- 
hicular activities (EVA’s) and evaluating EVA worksites 
is whether the EVA crew can reach and comfortably 
work at the designated worksite. Reach considerations 
are important because, in microgravity, astronauts must 
restrain their bodies in order to exert forces or moments 
on other objects. In the United States space program, 
body restraint during EVA is usually accomplished by 
foot restraints, particularly when tasks require applying 
significant forces. The articulated portable foot restraint 
(APFR) used on the international space station (ISS) can 
be adjusted to accommodate many orientations but can 
only be attached to the ISS’s structure at discrete loca- 
tions. Because of the constraints on foot restraint place- 
ment, it is essential to determine in the planning process 
whether a work site is reachable from an available foot 
restraint location. Although foot restraints are not the 
only available method of providing body restraint, the 
necessity of restraining the body for microgravity work 
makes the problem of providing both adequate restraint 
and access to the worksite a universal one in micrograv- 
ity, both inside and outside the spacecraft. A work enve- 
lope analysis that incorporates the mechanics of the 
space suit is a useful method for both assessing potential 
worksite locations and evaluating the functional signifi- 
cance of modifications to space suit mobility and visibil- 
ity. 

A reach envelope is the region in three-dimensional 
space that a person can reach. The work envelope is a 
subset of the reach envelope, representing the volume in 
which a person can comfortably work. Reach and work 
envelopes depend on the size and flexibility of the indi- 
vidual. A standard practice is to size workspaces to 
accommodate the reach and work envelopes of individu- 
als at the extremes of the expected size range, for exam- 
ple, 5th percentile females and 95th percentile males. 29 
Most work envelopes are determined experimentally, by 
measuring how far people of different sizes can reach 
and obtaining subjective information about the difficulty 
of working with the hands in different locations. The 
results of empirical work envelope studies are typically 
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a set of boundaries within which a specified percentage 
of the intended population can work comfortably. An 
alternative approach to reach envelope analysis has been 
developed recently, which uses robot kinematic analysis 
methods to determine the boundaries within which a 
person can reach with prescribed limits on joint ranges 
of motion. 15, 17, 36 Reach envelopes that are determined 
computationally can be tailored to specific individuals 
or populations of people. 

This report describes a novel work envelope analysis 
method, which combines inverse kinematics and a 
model of space suit mobility, enabling the work enve- 
lope of a space-suited astronaut to be assessed in greater 
detail and with more generality than experimental meth- 
ods allow. The computational work envelope analysis 
extends the inverse kinematics reach envelope work by 

prior investigators 15, 17, 36 to determine feasible working 
positions, which is a more restrictive category than 
reachable positions, incorporating human strength limits 


and visibility considerations. The model of space suit 
mobility developed in this work, which is derived from 
experimental data, predicts the torques needed to main- 
tain static positions of the space suit arm joints. 33 ' 34 

The work envelope analysis is based on previously 
described space suit joint mobility experimental data as 
well as incorporates the torque model. Torques required 
to maintain the space suit shoulder flexion, shoulder 
abduction, and elbow flexion joints were calculated 
using a linear regression model, which predicts torque 
as a function of joint angle and angular velocity, as 
shown in Equation 25. 

Torque = k\+ k 2 * angle + k 3 x angular velocity Eq 25 

Values of k], k 2 , and k 3 and the regression statistics R 2 
and standard error of the fitted torque are given in Table 

9 33 


Table 9: Torque model coefficients 33 


Joint 

k, (Nm) 

k 2 (Nm/ 
deg) 

k 3 (Nm/ 
(deg/sec)) 

R 2 

standard 
error (Nm) 

Shoulder 

flexion 

-23.4 

0.222 

0.120 

0.777 

4.39 

Shoulder 
abduction 
>50 deg 

-75.5 

1.20 

0.372 

0.743 

6.42 

Shoulder 
abduction 
<50 deg 

-27.2 

0.303 

0.110 

0.899 

0.553 

Elbow 

flexion 

-9.04 

0.155 

0.0380 

0.823 

2.12 


5.1 Work envelope criteria 

Three criteria detennine whether a proposed hand posi- 
tion is in the work envelope: the torques that are needed 
to hold the required arm configuration, visibility, and the 
shape of the resulting work envelope boundaries. The 
shape of the work envelope is considered because highly 
convoluted work envelope boundaries are difficult to 
use, since they are excessively sensitive to the relative 
positions of the EVA astronaut and the worksite. How- 
ever, adding the shape criterion overconstrains the prob- 
lem; as a result, a single torque limit cannot be enforced 
uniformly. To address this issue, a smoothing algorithm 


is incorporated to systematically trade off the joint 
torque and work envelope shape criteria. 

Joint torque limits 

In order for a person to work for several minutes in a 
particular arm configuration, the torques needed at all of 
the arm joints should be low enough that fatigue does 
not prevent the person from maintaining the desired arm 
position. The strength measure of interest in a work 
envelope analysis is the maximum torque that can be 
sustained for several minutes in a static position. In a 
survey of research on muscle fatigue, Astrand reported 
that the maximum force in static, isometric (constant 
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position) muscle contractions can be held for only a 
short period of time, but 1 5% of the maximum force can 
be held for at least 10 minutes, 30% of the maximum 
torque can be held for approximately 3 minutes, and 
50% of the maximum force can be held for 1.5 min- 
utes. 3 


“Low torque” and “excessive torque” limits are set for 
each joint and direction based on the relationship 
between maximum holding time and percentage of max- 
imum isometric torque. An arm position for which no 
joint torque exceeds 15% of the maximum torque is 
considered “low torque”. An arm position which 
requires more than 30% of the maximum torque from 
any joint is considered “excessive torque”. 

Visibility 

The space suit helmet and the display and control mod- 
ule (DCM), which is located on the front of the space 
suit, limit a space-suited astronaut’s field of view. A 
suited astronaut’s field of view is limited in the upward 
and left-right directions by the space suit helmet and 
visor and in the downward direction by the DCM and 
the mini-workstation, which is a tool rack located on the 
front of the hard upper torso (HUT). Only points that are 
within the field of view are included within the work 
envelope. 

Field of view limits for the EMU space suit published in 
a NASA EVA interface requirements document entitled 
NSTS-07700 1 were used in the work envelope analysis. 
Although the mini-workstation obstructs downward vis- 
ibility, the NSTS-07700 visibility limits do not include 
the presence of the mini-workstation. Consequently, the 
analysis presented herein does not consider the down- 


Isometric strength data for shoulder flexion, shoulder 
abduction, humerus rotation, and elbow flexion motions 
in both positive and negative directions, which are 
defined in Figure 3, are shown in Table 10 s . Torque lim- 
its used in the work envelope analysis are based on the 
limits in Table 10, 


ward visibility limitations imposed by the mini-worksta 
tion. 

70 deg 

f 

62 deg 



♦ 

70 deg 

FIGURE 34. Extravehicular Mobility Unit field of 
view 1 

Work envelope shape 

Using an arbitrary torque limit level to determine 
whether points are inside the workspace results in a 


TABLE 10. Isometric strength limits for shoulder and elbow. 5 


Motion 

50th 

percentile 
male (Nm) 

95th 

percentile 
male (Nm) 

50th 

percentile 
female (Nm) 

95th 

percentile 
female (Nm) 

Shoulder flexion (+) 

71 

101 

37 

57 

Shoulder flexion (-) 

67 

115 

30 

54 

Shoulder abduction (+) 

67 

103 

33 

57 

Shoulder abduction (-) 

92 

119 

40 

60 

Elbow flexion (+) 

77 

111 

41 

55 

Elbow flexion (-) 

46 

67 

27 

39 
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work envelope with a complex and convoluted shape. A 
work envelope with such a complicated shape is difficult 
to implement in planning EVAs or designing EVA work- 
stations. Small changes in the relative positions of the 
EVA astronaut and worksite may shift the worksite into 
or out of the work envelope. Excessive sensitivity to 
positioning is unwarranted because the appropriate max- 
imum torque levels depend on the strength and fatigue 
level of the astronaut, the duration of the task, and the 
accuracy of arm motions required in the task. Further- 
more, a complicated work envelope boundary requires 
storing a large number of boundary vertices for compu- 
tational implementation, and is difficult to visualize. 
Therefore, the boundaries of the work envelope should 
be as simple as possible, while including the maximum 
number of low-torque points and excluding all exces- 
sive-torque points. 

The shape of the work envelope is set by categorizing 
possible points, based on visibility and required torques, 
as “always include”, “always exclude”, and “possibly 
include”. Points that are visible and require torques that 
do not exceed a specified low torque limit are in the 
“always include” category, points that are not visible or 
require any joint torque to exceed a specified excessive 
torque limit are in the “always exclude” category, and 
points that are visible and require torques that fall 
between the low torque and excessive torque limits are 
in the “possibly include” category. The initial boundary 


is drawn to include only the “always include” points. 
Then the boundary is smoothed by adding a limited 
number of points that are in the “possibly include” cate- 
gory. Work envelope boundary smoothness is assessed 
by determining how many vertices are needed to draw a 
piecewise-linear boundary around the work envelope 
region using the contourc algorithm in Matlab (The 
Mathworks, Natick, MA). Minimizing the number of 
vertices needed to draw the boundary results in the 
smoothest boundary. 

Arm model 

The arm model chosen for the inverse kinematics calcu- 
lations is a simplified approximation of the complex 
mechanics of the human arm. The number and type of 
articulations of the model arm were chosen to coincide 
with the joints included in the space suit model and to 
limit the kinematic redundancy of the model arm. The 
model arm includes four joints: shoulder flexion (sf), 
shoulder abduction (sa), humerus rotation (hr), and 
elbow flexion (ef). 

The lengths of the arm segments are set based on 
anthropometric data 29 , which is listed in Table 1 1 . The 
lower arm segment length is the elbow to wrist distance 
plus half the hand length. Arm segment lengths can also 
be customized to a specific person. 


TABLE 11. Arm segment lengths 29 


Segment 

Male, 50th 

Male, 95th 

Female, 50th 

Female, 95th 

percentile (cm) 

percentile (cm) 

percentile (cm) 

percentile (cm) 

Upper arm 

36.6 

39.4 

29.8 

32.4 

Lower arm 

38.3 

41.6 

33.1 

35.3 


5.2 Work envelope methods 

The visibility and joint torque limitations are applied in 
a five-step process. First, visibility is evaluated: if the 
target point is outside of the field of view, it is excluded 
from the work envelope. Second, inverse kinematics is 
used to determine the joint angles at the shoulder and 
elbow that are needed to place the hand on target. The 
inverse kinematics solution is not unique, because three 
coordinates of hand position determine four joint angles. 
The inverse kinematics analysis outputs a set of several 
arm configurations that put the hand on target. Third, the 
space suit torque-angle model accepts these joint angles 
as inputs and calculates the joint torques required to 
maintain the specified arm configuration. Fourth, the 


joint torques are compared to the pre-set torque limits 
and the arm configurations for which no joint torques 
violate the torque limits are returned as feasible arm 
positions. If at least one feasible arm position exists for 
a target hand position, the target point is considered to 
be within the work envelope. Finally, the overall diffi- 
culty of reaching the target point is assessed by calculat- 
ing a reach difficulty metric, based on the required joint 
torques. The feasible arm configuration that has the low- 
est reach difficulty metric is recorded for each target 
position in the work envelope. 

Visibility 

Visibility of the target point is evaluated first, because 
the visibility constraint eliminates a large portion of the 
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possible target points. Visibility is evaluated by locating 
the suited person’s eyepoint relative to the right shoul- 
der, then using the eyepoint position and target hand 
position to determine whether the target hand position 
lies within the angular field of view limits defined in 
Figure 34. A single point that represents the eyepoint 
position is located on the body centerline, halfway 


TABLE 12. Eyepoint position relative to the right shoulder 29 


Percentile 

M 50 th 

M 95th 

F 50 th 

F 95th 

Vertical 

(cm) 

16.5 

16.9 

15.9 

16.5 

Left-right 

(cm) 

20.6 

22.2 

18.7 

20.3 


between the eyes. The eyepoint coordinates that were 
used for 50th and 95th percentile males and females are 
listed in Table 12. No anthropometric data was avail- 
able for the position of the eyes forward of the shoulder, 
therefore, a distance of 10 cm was used in all cases. 


Inverse kinematics 

Inverse kinematics refers to the process of determining 
the joint angles that are required to place an end effector 
at a specified point in space. Inverse kinematic analysis 
of the arm is calculated using a geometric method 

described by Korein 19 for a four degree of freedom arm. 
The inverse kinematics analysis was carried out for a 
three-dimensional mesh of points, evenly spaced 6.2 cm 
apart, representing proposed hand positions. 


The 4 degree of freedom arm is shown in Figure 35. 
Points S, E, and W represent the shoulder, elbow, and 
wrist joint locations. Coordinates are centered on point 
S, which is at the center of the right shoulder, and coor- 
dinate directions are defined in relation to the torso. The 
z axis points down, the y axis points forward, and the x 
axis points to the left. 



FIGURE 35. A.Two-link arm, showing the shoulder 
(S), elbow (E) and wrist (W) locations. 
B. Planar view of two-link arm. 

The hand target position uniquely determines the elbow 
angle, while coordinated displacements of the shoulder 
flexion, shoulder abduction, and humerus rotation joints 
allow the arm to assume multiple configurations that all 
place the hand on the target. If 360 deg rotation of all 
joints were allowed, the locus of possible elbow posi- 


tions that place the hand on target would be the dashed 
circle shown in Figure 35A. Since the joint ranges are 
limited, the elbow position locus forms arcs on the cir- 
cle, like the solid arc in Figure 35A. In virtually all 
cases, there are many possible arm configurations that 
place the hand at its target position. 

The following procedure, adapted from Korein’s 19 
human arm inverse kinematics method, is used to deter- 
mine if a proposed hand target point W lies within the 
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arm reach space and calculate the shoulder flexion, 
shoulder abduction, humerus rotation and elbow flexion 
angles required to place the hand at the target position. 
Rot(x,0) indicates the rotation matrix required to rotate 
a vector about an axis described by vector x by angle 

ft 15 

1 . Determine the elbow angle e. 


The elbow angle e is determined uniquely by the dis- 
tance between the shoulder joint and the hand target 

point, |wj . Based on |w| and the lengths of the upper 
and lower arm L u and Lj, e is calculated from the law of 
cosines, as shown in Figure B and Equation 26. 


e = jt - acos 


'( L l + L l - 


M 2 )! 


2L..L 


u^L 


Eq 26 


The hand position W| that results from the elbow joint 
rotation alone is determined by 27. The hand position 
Wj is a three-element vector in the coordinate frame 
given in Figure 35A. 



0 


sin(jt - e) 

*1 = L u 

0 

+ l l 

0 


_ 1 _ 


-cos(ji - e) 


2. Prescribe a humerus rotation angle, t. 


Eq 27 


from w 2 and w, in order to rotate w 2 to w. The unit vec- 
tor in the x-y plane that is equidistant from w 2 and w is 

defined as a and is given by: 


-( w x~ w x2 ) 


0 


( w ’.v _w ’y 2)1 

K- w *2) 

0 


Eq 29 


The shoulder rotation angle required, ?/, is the angle 
between the projections of w 2 and w on the plane per- 
pendicular to a . To determine q, the projections of w 2 

and w on the plane perpendicular to a are defined as 
follows. 


pw = (ax w) x a 
pw , = (ax w 2 ) x a 


Eq 30 


The angle between pw and pw 2 determines q. 


q = acos 


pw • pw 2 
|pw||pw 2 


Eq 31 


A range of humerus rotation angles, r, distributed 
evenly over the available humerus rotation range, -113 
deg to 70 deg, is chosen. The hand position that results 
from each of the humerus rotation angles is then calcu- 
lated. 


w 2 = Rot(z,T)vV| Eq 28 

Any choice of rthat results in permissible values for the 
other three arm joint angles (shoulder flexion, shoulder 
abduction, and humerus rotation) is an allowed arm con- 
figuration. To resolve the kinematic redundancy, the 
inverse kinematics calculations were carried out in par- 
allel for multiple values of r, then one arm configuration 
was chosen to minimize the torques required at the 
joints to maintain the arm configuration. 

3. Rotate the shoulder to bring w 2 to w, without chang- 
ing the humerus rotation angle. 

The shoulder flexion and shoulder abduction rotations 
are performed simultaneously by rotating vector w 2 
about the shoulder so that the hand moves to its target 
location. In order to preserve the r specified in step 2, 
the rotation axis must be chosen to avoid rotations about 
the z axis. The rotation axis must also be equidistant 


Since the numerator and denominator in Equation 3 1 are 
both always positive, the sign of q is determined accord- 
ing to Equation 32, so that the rotation is right-handed 

with a . 


q > 0, w x m’ 2 • a > 0 
»1 < 0, w x w 2 * a < 0 


Eq 32 


This single rotation of angle i] about axis a will later be 
decomposed into shoulder flexion and shoulder abduc- 
tion angles. 

4. Determine the elbow position. 

The new elbow position g| that results from the shoulder 
rotation is determined by rotating the vector from point 
S to point E by angle q about axis a . 

e { = eRot(§, q) Eq 33 

5. Determine the shoulder flexion and shoulder abduc- 
tion angles. 

The shoulder flexion and shoulder abduction angles can 
be determined by the vector from the shoulder S to the 
elbow £[ . The shoulder flexion and shoulder abduction 
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degrees of freedom are assumed to be successive rota- 
tions about perpendicular axes. The shoulder flexion 
(SF) rotation, which is performed first, is a rotation 
about the negative x axis. The shoulder abduction (SA) 
rotation, performed second, is a rotation about the a neg- 
ative y axis that is attached to the upper arm segment 
and rotated with the shoulder flexion rotation. While the 
order of the two shoulder rotations must be specified 
and applied uniformly, there is no inherently preferred 
ordering. 



At this point in the procedure, angles for the four arm 
joints have been calculated for ten arm configurations 
corresponding to the ten prescribed values of humerus 
rotation angle. All of these arm configurations place the 
hand on the target, although some of them are unattain- 
able because they exceed the available joint ranges of 
motion. The space suit model is used next to determine 
the torque required to maintain each of these arm con- 
figurations. A single arm configuration is then chosen as 
a unique solution, based on the torques needed to hold 
the configuration. 

Space suit joint torques 

Torques required to hold the space suit in the configura- 
tions determined by the inverse kinematics analysis 
were calculated from the regression model given in 
Equation 25 and Table 9. Because the arm position is 
considered to be static in the work envelope configura- 
tion, the angular velocity term in Equation 25 is set to 
zero, to isolate the torque needed to hold a specific joint 
angle. 

Reach difficulty metric 

The difficulty of reaching a specific location is indicated 
by the reach difficulty metric, which is a nondimen- 
sional, positive scalar corresponding to each proposed 
arm configuration. A low reach difficulty metric means 
that the arm configuration is easy to maintain; a high 
reach difficulty metric means that the arm configuration 
is difficult to maintain. The reach difficulty metric is the 
sum over three arm joints: shoulder flexion, shoulder 
abduction, and elbow flexion, of the required joint 
torque divided by either the positive or negative torque 
limit for that joint, as shown in Equation 35. The 
humerus rotation joint is not included because torque is 


not required to hold a static humerus rotation joint posi- 
tion. 


joint 3 

metric = V ( — torque \ Eq 35 

lj \ torque limit/ 

joint 1 

For each proposed hand location, the inverse kinematics 
analysis outputs several arm configurations that place 
the hand at the target location. The torques required to 
maintain these arm configurations are compared to the 
torque limits. If the torque at any joint exceeds its limit, 
that arm configuration is eliminated from consideration. 
When more than one ami configuration is allowable, the 
one with the lowest reach difficulty metric is chosen as 
the single solution for that hand location. Because a sin- 
gle arm configuration is chosen for each proposed hand 
location, one value of the reach difficulty metric corre- 
sponds to each point in the proposed workspace. There 
are no constraints applied to the reach difficulty metric; 
the torque limits that detennine whether a point is 
within the work envelope apply only to torques required 
from individual joints. 

Boundary smoothing 

The boundaries of the work envelope regions are 
smoothed using image processing techniques to set the 
simplest work envelope boundaries that include as many 
low-torque points as possible and exclude excessive- 
torque and nonvisible points. The smoothing technique 
chosen is objective and systematic and uses standard 
Matlab library functions, although other smoothing 
techniques may be used to accomplish similar results. A 
binary image is formed from a horizontal or vertical pla- 
nar section of the workspace, with points in the “always 
include” category set to be black and points outside that 
category set to be white. The “always include” category 
is the region of points that is visible and requires no 
more than 15% of the maximum torque for each joint. 
Additional points from the “possibly include” category 
are added to this region, or colored black, using mor- 
phological operators to fill interior gaps in the black 
region. A gaussian low-pass filter is then applied to the 
image, resulting in a blurred, grayscale version of the 
image. The grayscale image is converted back to binary 
form, using a variable threshold value. The threshold is 
chosen to meet two conditions: 

1 . The smoothing process may add no more than 30% 
new points to the image 

2. A piecewise linear boundary drawn around the white 
region in the smoothed image should require the 
minimum number of vertices. 


Quantifying Astronaut Tasks: Robotic Technology and Future Space Suit Design - PI: Dava Newman 


35 



The smoothing algorithm was carried out first for hori- 
zontal slices, then for vertical slices through the work- 
space region. 


5.3 Results 

Smoothed work envelope boundaries Three- 
dimensional plots of the work envelope boundaries that 
were calculated using female and male anthropometric 
data for 50th percentile size and 95th percentile strength 
are shown in Figure 36 and Figure 37. 



60 



FIGURE 36. Smoothed work envelope boundaries 


for female, 50th percentile size, 95th 
percentile strength. 
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FIGURE 37. Smoothed work envelope boundaries 
for male, 50th percentile size, 95th 
percentile strength. 



Horizontal slices through the work area, illustrating the 
reach difficulty metric and the smoothed work envelope 
boundaries are shown in Figure 38 for the female 50th 
percentile size, 95th percentile strength case and 
Figure 39, for the male, 50th percentile size, 95th per- 
centile strength case. 
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FIGURE 38. Reach difficulty metric shown by 

intensity and work envelope boundaries 
shown by black line for female, 50th 
percentile size and 95th percentile 
strength. A.) 6 cm below shoulder. B.) At 
shoulder level. C.) 6 cm above shoulder. 
D.) 19 cm above shoulder. 
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FIGURE 39. Reach difficulty metric shown by 

intensity and work envelope boundaries 
shown by black line for male, 50th 
percentile size and 95th percentile 
strength.A.) 37 cm below shoulder. B). 
19 cm below shoulder. C.) At shoulder 
level. D.)19 cm above shoulder. 

Joint angle ranges The space suit torque model, 
arm segment lengths, and joint torque limits determine 
the angle range that is generated by the work envelope 
analysis for each joint. Ranges of motion of each of the 
four arm joints for torque limits of 15%, 30%, 50% and 
1 00% of the maximum isometric strength are shown in 
Figure and Figure . These ranges of motion are the 
maximum and minimum joint angles selected in the 
inverse kinematics analysis. While they do span the 
majority of available joint angles, they are not mathe- 
matically guaranteed to cover the entire range of feasi- 
ble joint angles predicted by the space suit torque model 
and the strength limits. Consequently, small numerical 
differences can be seen between the angle ranges for 
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cases that have the same strength limits but different arm 
lengths. 


200 



sh flex sh abd 

FIGURE 40. Joint angle ranges for female, 50th 

percentile size, 95th percentile strength 
deduced from work envelope analysis. 
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FIGURE 41. Joint angle ranges for male, 50th 

percentile size, 95th percentile strength 
deduced from work envelope analysis. 

Work envelope volume The volume enclosed 
within the work envelope boundaries serves as one glo- 


hum rot elb flex 

bal measure of work envelope size, for purposes of com- 
parison. The cases tested in the work envelope analysis 
produced different volumes for the boundaries enclosing 
the 15%, 30%, 50% and 100% strength contours. Vol- 
umes enclosed by these boundaries are shown in 
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Figure 42. Work envelope volume is shown in 
Figure 43, nonnalized by the volume of a hemisphere 
whose radius is equal to the arm length for the particular 
case tested. 
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FIGURE 42. Volume enclosed within the 15%, 30%, 
50%, and 100% torque limit boundaries 
for different anthropometric conditions. 
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FIGURE 43. Volume enclosed within the 1 5%, 30%, 
50%, and 100% torque limit boundaries 
for different anthropometric conditions, 
nonnalized by available volume. 

Sensitivity analysis 

The three work envelope criteria, torque limits, visibil- 
ity, and boundary shape, contribute to the size and shape 
of the work envelope. Because multiple constraints set 
the boundaries of the work envelope, it is informative to 
determine which constraints affect the outcome of the 
work envelope analysis and which constraints do not. 
The combination of space suit joint stiffness and human 
strength limits sets the range of motion for each joint. 
The joint ranges of motion, combined with the lengths 
of the arm segments, determine the region that a person 
can reach. In addition, visibility limits may exclude 
reachable but out of view areas from the work envelope. 
Using anthropometric data that represents the size and 
strength of real people makes the relationship between 
the work envelope criteria and the resulting work enve- 
lope more complicated. For instance, increasing a high, 
but realistic, torque limit may not increase the size of the 
work envelope if the joint can already traverse its full 
range with the lower torque limit. The results of the sen- 
sitivity analysis illustrate how the astronaut’s size, 
strength and visibility and the space suit’s stiffness 
affect the size and shape of the work envelope. 


Human body size Body size enters the inverse 
kinematics analysis through the lengths of the upper and 
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lower arm segments. The joint angle ranges of motion 
shown in Figure 40-Figure 41 indicate that body size 
had only a very small effect on the ranges of motion of 
the joints. The small changes that are seen in joint 
ranges of motion between 50th percentile size and 95th 
percentile size people occur because the inverse kine- 
matics analysis does not guarantee use of all of the 
available range of the joints. Because changing the arm 
segment lengths changes the required joint angles for a 
given hand position, it is expected that the joint angles 
output from the inverse kinematics analysis will be 
slightly different for different arm lengths. 

The work envelope volume, plotted in Figure 42, and 
the percentage of available volume occupied by the 
work envelope, shown in Figure 43, give an overall pic- 
ture of the work envelope size. At the 50% and 1 00% 
torque limits for the male anthropometric data and the 
female 95th percentile strength anthropometric data, the 
work envelope volume is determined solely by the arm 
length. The dominant effect of arm length on work enve- 
lope volume at the female 95th percentile strength and 
male 50% and 100% torque limits can be seen in 
Figure 43, where the percentage of available volume 
occupied by the work envelope is virtually the same in 
all cases. The uniformity of work envelope volumes for 
the male, 50% and 100% torque limits is explained by 
Figure 41, which indicate that the joint ranges of motion 
for 50% and 100% torque limits vary only slightly over 
all four male size and strength combinations. A different 
trend is seen for the female, 50th percentile strength, 
anthropometric data, where work envelope volume con- 
sistently increases as torque limits increase. 

When the torque limits are high enough that the full 
ranges of the joints are available, the arm length deter- 
mines the work envelope volume. This is the case for the 
50% and 100% male torque limits and the 50% and 
100% torque limits for female, 95th percentile strength, 
but it is not true for the female 50th percentile strength 
data, because the female 50th percentile strength torque 
limits are lower. Both arm length and strength limits 
affect the work envelope size for the female 50th per- 
centile strength anthropometric data, at all torque limit 
levels. 

Human strength limits and space suit joint 
stiffness Human strength limits and the torques 
required to bend the space suit joints together determine 
the available range of motion of each joint. As the 
strength limits increase or the space suit joints requires 
less torque, the range of motion expands. Benefits asso- 
ciated with increasing torque limits or decreasing space 
suit stiffness are exhausted when the joints achieve their 
full range. When the full range of a joint is available, 


increasing torque limits further does not increase the 
size of the work envelope, and the work envelope size 
and shape are limited by the visibility constraint and 
arm segment lengths. Changes in work envelope volume 
due to changes in torque limit levels ( 1 5%, 30%, 50% or 
100%) and strength percentile values (50th percentile or 
95th percentile) illustrate the effects of across-the-board 
reductions in space suit stiffness, while a more detailed 
analysis of the effects of increasing individual joint 
torque limits shows the effects of improving mobility in 
single joints. 

Figure 43 shows the normalized volume occupied by the 
work envelope for 15%, 30%, 50% and 100% torque 
limits. For the male, 50th percentile size, 95th percentile 
strength case, increasing the torque limits from 1 5% to 
30% results in an increase in work envelope volume of 
270%, while increasing the torque limits further to 50% 
results in an additional gain in work envelope volume of 
only 5%. Therefore, reducing the stiffness of all space 
suit joints by a factor of two results in a large increase in 
work envelope volume, while further reductions in 
space suit stiffness produce little gain. 

Looking at the effects of increasing torque limits for 
individual joints indicates which joints limit the overall 
size of the work envelope and where each joint's limits 
bound the work envelope. Individual joints were ana- 
lyzed using the female, 50th percentile size, 95th per- 
centile strength case. The baseline case uses torque 
limits of 30% of the maximum strength for all joints. Six 
variations on the baseline case were analyzed, varying 
the positive and negative torque limits on the shoulder 
flexion, shoulder abduction, and elbow flexion joints. 
The humerus rotation joint was not considered because 
the space suit’s arm bearing allows it to traverse its full 
unsuited range of motion without reaching the 30% 
torque limits. In each variation evaluated in the sensitiv- 
ity analysis, the positive or negative torque limit for one 
joint was increased from 30% of the strength limit to 
50% of the strength limit. The work envelope volume in 
the baseline case and the six variations, normalized by 
the volume of a hemisphere with radius equal to the total 
arm length, are shown in Figure 44. The largest gain in 
work envelope volume is associated with an increase in 
the negative shoulder abduction torque limit. Increasing 
the negative shoulder flexion torque limit also increases 
the work envelope volume somewhat, while positive 
shoulder flexion and positive and negative elbow flexion 
do not limit the work envelope volume. 
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FIGURE 44. Percentage of available volume 

occupied by work envelope when torque 
limits are increased from 30% to 50% 
for single joints. The baseline case uses 
torque limits of 30% of the strength 
limits for all joints. 

Visibility The visibility constraint excludes areas 
close to the body, at the top and bottom of the work 
space, from the work envelope. The 30% torque con- 
tours calculated using field of view restrictions were 
compared to boundaries calculated using the same 
torque limits but no field of view restrictions for male 
50th percentile size, 95th percentile strength and female 
50th percentile size, 95th percentile strength. As 
Figure 45 shows, the visibility limits exclude areas from 
the work envelope at the top and bottom, close to the 
body, but they do not change the right and left sides of 
the work envelope boundaries. Removing the visibility 
constraint increases the volume occupied by the work 
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envelope by 33% in the male case studied and 22% in 
the female case. 
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FIGURE 45. Comparison of 30% torque contours 
with and without visibility limits, for 
female, 50th percentile size, 95th 
percentile strength. A.) 12 cm below 
shoulder. B.) 19 cm above shoulder. C.) 

3 1 cm above shoulder. D.) 44 cm above 
shoulder. 

The work envelope analysis demonstrates the practical 
utility of modeling the space suit mechanics on a funda- 
mental level, yielding both operational and design rec- 
ommendations. Using models of the torque-angle 
relations of the space suit joints and space suit kinemat- 
ics, data on visibility in the space suit, and a spatial 
smoothing algorithm, a rapid, systematic, and customiz- 
able method was developed to determine work envelope 
boundaries for a space-suited astronaut.The method was 
used to calculate work envelope boundaries for two 
cases, using both male and female anthropometric data. 
In addition to boundaries on suitable locations of work 
areas, the work envelope analysis also indicates the rela- 
tive difficulty of holding the hand in different locations 
while wearing a space suits. 


Two important simplifying assumptions were made in 
the analysis, which limit the numerical and statistical 
accuracy of the results. First, the joints of the human 
arm were assumed to have constant rotation axes 
throughout their entire range of motion, resulting in 
potential errors in computed hand position of approxi- 
mately 1 cm-2cm. These hand position errors are within 
the 6.2 cm mesh spacing for trial hand positions. Sec- 
ond, several human body dimensions, which were 
obtained from anthropometric databases at the 50th per- 
centile and 95th percentile levels, were combined in this 
analysis. Because the anthropometric data were com- 
bined, the work envelopes generated by the analysis 
should not be considered to be exact 50th percentile or 
95th percentile work envelopes. 

A sensitivity analysis of the work envelope boundaries, 
which varied joint torque limits, arm length, and visibil- 
ity limits, yielded several space suit design recommen- 
dations. Increasing mobility of the space suit’s shoulder 
joint to allow the arm to be brought closer to the body 
centerline in front of the body results in the largest gains 
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in work envelope size. An across-the board reduction in 
space suit joint stiffness by a factor of 2 results in large 
gains in work envelope size, while further reductions in 
space suit joint stiffness only yield small gains. Visibil- 
ity improvements in the upward and downward direc- 
tions increase the work envelope size, while visibility 
does not limit work envelope size in the left and right 
directions. 

6.0 Summary and Future EVA 
Research 

The primary aim of this report was to summarize the 
research effort undertaken to advance the current under- 
standing of astronauts’ capabilities and limitations in 
space-suited EVA by compiling a detailed database of 
the torques needed to bend the joints of a space suit, 
developing models of the mechanics of space suit joints 
based on experimental data, and utilizing these models 
to estimate a human factors performance metric, the 
work envelope for space suited EVA work. 

6.1 Space suit mobility database 

The space suit torque-angle database was collected 
using an Extravehicular Mobility Unit space suit, in a 
novel measurement approach which used both human 
subjects and an instrumented robot to collect joint angle 
and joint torque data for realistic human motions. 
Human test subjects carried out arm and leg motions 


both wearing the space suit and not wearing the space 
suit, supplying realistic joint angle trajectories for each 
of 20 motions. The joint angle trajectories produced by 
the human subjects were then used as command inputs 
for the instrumented robot, so that the robot imitated the 
humans’ motions while torques were measured at each 
of the robot’s joints. Torques on the robot’s joints due to 
the weight of the robot and space suit were subtracted 
from the torques that the robot measured, resulting in a 
consistent set of joint angle and torque data. 

The database compiled in this work is more extensive 
than any other published space suit torque-angle data 
set, covering 1 1 joints over a large range of angles, 
including multiple degrees of freedom at the shoulder, 
hip, and ankle. Realistic, three-dimensional human-gen- 
erated motions were used for data collection. Because 
torque data was collected by the robot as a surrogate for 
a human occupant of the space suit, the torques mea- 
sured in this study are more representative of realistic 
conditions than data from previous studies 2, 7 ’ 27 that 
measured joint stiffnesses for empty, pressurized space 
suits. A limited comparison of torques required to bend 
the elbow and knee to two angles, shown in Table 13, 
illustrates that the torques measured in this study are 
higher than the torques measured with empty space 
suits. The space suit torque-angle database serves as a 
basis for developing and validating both mathematical 
and physical models of space suit joint mobility charac- 
teristics. 


TABLE 13. Space suit joint torques 


Study 

Dionne 

Menendez et 
al 

Abramov et al. 

Morgan et al. 

Current Study 

Methods 

EMU 

empty 

suit 

Prototype 
joint seg- 
ments, 
empty 

Orlan-DMA, 
4.3 psi 
(30 KPa) 
empty suit 

EMU 

human 

subjects 

EMU, human 
subjects and 
robot 

Knee, 72 deg. 

3.2 Nm 

NA 

6.0 Nm 

8.1 Nm 

3.74±.0676 Nm 

Elbow, 80 deg 

2.0 Nm 

2 Nm 

2.2 Nm 

3.4 Nm 

14.6±.136 Nm 


6.2 Modeling 

The space suit torque-angle database presents a rare 
opportunity to develop models of space suit mobility 
and verify them against experimental data. Models of 
space suit mobility are useful in two applications: first, 
numerically predicting the torque required to bend a 
space suit’s joints and second, understanding the physi- 


cal processes that determine the stiffness of a space 
suit’s joints. These two applications clearly require two 
different modeling approaches: a descriptive mathemati- 
cal modeling technique based on experimental data and 
a theoretical model based on physical principles. Sec- 
tion 4 describes two modeling efforts: a mathematical 
model based on empirical data that predicts the torque 
required to bend the space suit’s joints and a comparison 
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between two physics-based models of bending pressur- 
ized cylinders and the space suit torque angle database. 

Numerical predictions of the torque required to bend 
space suit joints were generated using a mathematical 
hysteresis modeling technique. The experimental data 
reported in Section 3, as well as other investigators’ 
data, indicate that the torque-angle relationship for 
space suit joints is markedly hysteretic. The Preisach 
hysteresis model, which has been used to model other 
hysteretic mechanical systems, is used to mathemati- 
cally model the torque required to bend the space suit 
joints as a function of bending angle history. The Prei- 
sach model coefficients were identified for the torque- 
angle characteristic of six of the space suit’s joints: 
elbow flexion, hip flexion, hip abduction, knee flexion, 
ankle rotation and ankle flexion. In addition, a new 
method was developed for estimating the variance of the 
error in the Preisach model’s torque predictions. Com- 
parisons between the Preisach model’s output and 
experimental data indicate that the Preisach hysteresis 
model accurately reproduces torques needed to bend 
space suit joints. 

The physics-based modeling effort compares two 
approximate models of pressurized fabric cylinder bend- 
ing to experimental data. The two models, the beam 
model and the membrane model represent limiting cases 
on the physical processes that govern space suit joint 
mobility. The beam model assumes that space suit joint 
bending occurs through elastic deformations of the 
space suit fabric, with no change in internal volume 
when joints are bent. In contrast, the membrane model 
assumes that the space suit fabric never stretches, and 
the work done to bend the space suit joint goes entirely 
into compressing the gas inside the space suit. Compar- 
ing these two approximate models to experimental data 
illustrates the relative importance of elasticity and vol- 
ume loss on space suit joint mechanics and indicates 
whether the space suit behavior can be approximated by 
one model or the other. Understanding the physical pro- 
cesses that govern space suit joint performance leads to 
insights into how space suit joints would perform if their 
construction or operating conditions were changed. 

Comparisons between the two models and experimental 
data show that the membrane model matches the elbow 
and knee joint torque-angle data well over a bending 
angle range 30-50 degrees from the equilibrium bending 
angle and underestimates the torque required to bend the 
space suit joints for more extreme bending angles. The 
beam model is inconsistent with both the magnitude and 
curve shape of the elbow and knee joint torque-angle 
data. Data from the hip abduction joint, which has no 
mobility features and is consequently very stiff, is con- 


sistent with the beam model, although the limited angle 
range of the data prevents the beam model from show- 
ing important qualitative characteristics. 

Because the torque-angle performance of the space suit 
joints can be accounted for by changes in internal vol- 
ume without stretching of the space suit fabric, the 
mobility of the EMU elbow and knee joints is not sensi- 
tive to changes in the material properties of the fabric. 
The membrane model predicts that the design parame- 
ters that affect space suit mobility are the internal pres- 
sure, radius of the space suit segment, and the pleat 
design that allows the joint to be bent without reducing 
the internal volume. 

The physics-based modeling work in this thesis provides 
insights into the physical processes that govern space 
suit mobility, points out promising methods of improv- 
ing space suit mobility and provides a theoretical basis 
for further experimentation to develop more sophisti- 
cated physics-based models of space suit mobility. 

6.3 The work envelope: Applying space suit 
modeling to EVA operations 

The work envelope analysis demonstrates the practical 
utility of modeling the space suit mechanics on a funda- 
mental level, yielding both operational and design rec- 
ommendations. Using models of the torque-angle 
relations of the space suit joints and space suit kinemat- 
ics, data on visibility in the space suit, and a spatial 
smoothing algorithm, a rapid, systematic, and customiz- 
able method was developed to determine work envelope 
boundaries for a space-suited astronaut.The method was 
used to calculate work envelope boundaries for two 
cases, using both male and female anthropometric data. 
In addition to boundaries on suitable locations of work 
areas, the work envelope analysis also indicates which 
areas are best for workstation placement. 

A sensitivity analysis of the work envelope boundaries, 
which varied joint torque limits, arm length, and visibil- 
ity limits, yielded several space suit design recommen- 
dations. Increasing mobility of the space suit’s shoulder 
joint to allow the arm to be brought closer to the body 
centerline in front of the body results in the largest gains 
in work envelope size. An across-the board reduction in 
space suit joint stiffness by a factor of 2 results in large 
gains in work envelope size, while further reductions in 
space suit joint stiffness only yield small gains. Visibil- 
ity improvements in the upward and downward direc- 
tions increase the work envelope size, while visibility 
does not limit work envelope size in the left and right 
directions. 
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6.4 Recommendations for future work 

Experimental work Joints in the space suit include 
both bending joints and rotational bearing joints, to 
allow motion in multiple degrees of freedom in the arm. 
While this thesis focuses mostly on the performance of 
the bending joints, experimental work to quantify the 
motion of the space suit shoulder segment between the 
two bearings in the arm, as distinguished from the joint 
motion of the person wearing the space suit, would lead 
to a better understanding of the multiple degree of free- 
dom motion at the shoulder joint. 

Future experimental work would also contribute to the 
physics-based modeling efforts. The membrane model, 
which agrees closely with experimental data for the 
elbow and knee, predicts that joint stiffness increases 
linearly with space suit pressure. This result could easily 
be tested by measuring joint torque vs. angle at a range 
of pressures. Other predictions that the membrane 
model makes, including inextensibility, the amount of 
internal volume change, and the bent shape of the space 
suit joint could be tested in experiments with pressur- 
ized mock-ups of individual space suit joints. 

The work envelope analysis makes detailed predictions 
of easily-reachable locations based on body size and 
strength. For these predictions to have practical value, 
the work envelope analysis should be experimentally 
validated. Experimental work to validate the work enve- 
lope predictions should use arm length, viewpoint loca- 
tion, and static joint strength data from the test subjects 
to predict work envelopes tailored to those individuals. 
Validation of the calculated work envelope should be 
performed in neutral buoyancy simulation, rather than a 
1-g lab environment, so that the test subjects are not 
required to supply additional joint torques to support the 
weight of the space suit’s arm. Test subjects would 
attempt to hold one hand at pre-determined locations for 
several minutes and subjectively rate the difficulty of 
holding the hand position, as well as indicating whether 
the point is visible. Testing points should be chosen to 
fall both inside and outside of the calculated work enve- 
lope, with particular attention paid to points that are 
inside the calculated work envelope but outside of the 
NASA work envelope. Data collected should include 
whether the subject was successful in holding the hand 
at the test point for the required duration, the length of 
time that the test subject was able to maintain the 
required hand position, the subjective difficulty ratings 
and whether the point was visible. This information 
should then be compared to the calculated work enve- 
lope’s boundaries, the reach difficulty metric, and the 
visibility predictions. Modifications to the computa- 
tional work envelope based on these evaluations may 


include changes to the percentages of maximum torque 
that are used as torque limits and more detailed visibility 
information. 

Analysis Future analytical work should include 
implementation of existing space suit models and devel- 
opment of new models, based on experimental data. We 
have previously demonstrated in simulations that 
torques induced by the space suit have an important 
impact on body positions and loads when astronauts 
move large objects in EVA 30 , 31 . Future EVA dynamic 
simulation studies should include the torques needed to 
bend the space suit joints established under this research 
effort. 

New, physics-based models of space suit joint mobility 
should address joint design issues by including details 
of space suit joint construction. Work should also be 
done to asses the proper way of relaxing the inextensi- 
bility constraint of the membrane model to include the 
limited fabric stretching that likely occurs when space 
suit joints are bent. Hysteresis in fabrics is generally 
thought to be caused by the fabric threads sliding over 
each other as the fabric stretches 13 . Although the space 
suit joints show significant hysteresis, implying that 
there is some internal relative motion in the fabric, the 
membrane model’s agreement with experimental data 
implies the contrary, that the space suit fabric is inexten- 
sible. It would be particularly interesting to account for 
the hysteresis mechanism in light of this apparent con- 
tradiction. 

EVA standards and requirements The NSTS- 

07700 EVA requirements document 1 specifies NASA’s 
requirements for the space suited work envelope and the 
field of view for an EVA astronaut (Figure 34). The cur- 
rent requirements are simply stated, requiring 4 numbers 
for the work envelope (3 coordinates of cylinder center 
location and cylinder diameter) and 5 numbers for the 
field of view (angular field of view in 5 directions). The 
simple requirements are easy to understand and visual- 
ize, lending themselves to qualitative evaluation and 
application to worksite design. However, much more 
powerful methods, including CAD modeling, are cur- 
rently available for EVA worksite analysis and have the 
capability to enforce complicated requirements quanti- 

i 1 1 

tatively. ’ Now that advanced techniques exist for 
evaluating work envelope and visibility requirements, 
the work envelope and visibility standards should be 
updated to a higher accuracy and level of detail. 

The work envelope standard in NSTS-07700 is signifi- 
cantly smaller than the work envelope calculated based 
on space suit mobility and human strength limits. In par- 
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ticular, there are large areas beyond the left and right 
boundaries of the NSTS-07700 work envelope that are 
inside the calculated work envelope. Also, a visibility 
analysis of the NSTS-07700 work envelope using the 
visibility requirements from the same standards docu- 
ment indicates that a portion of the NSTS-0770 work 
envelope is not visible. These disagreements should spur 
a re-evaluation of the NSTS-07700 work envelope, 
including visibility, body size differences, and, at the 
least, enlarging its left-right dimension. 

The visibility standard in NSTS-07700 specifies the 
angular field of view in 5 directions. The lack of spatial 
resolution prevents this guideline from accurately speci- 
fying the visibility limitations caused by the shape of the 
helmet and visors and the geometry of other obstruc- 
tions, including the space suit’s Display and Control 
Module (DCM), and the mini-workstation, which both 
block downward visibility. Since powerful analytical 
techniques currently exist for assessing adherence to 
these requirements, the visibility requirement should be 
updated to improve its accuracy and make full use of the 
available analytical techniques for assessing compliance 
to requirements. 

6.5 Applications to current EVA 

The analytical techniques and space suit mobility data- 
base developed in this thesis can be applied to current 
EVA planning, training and hardware design and future 
planetary EVA system design. The computational work 
envelope method can be used for EVA and neutral buoy- 
ancy simulation planning and EVA hardware design 
evaluations. A larger and more realistic work envelope 
estimate would make more efficient use of time and 
resources in EVA work than the current work envelope 
estimate. The mathematical model developed in this the- 
sis is implementable in real time and suitable for incor- 
poration into EVA dynamic simulations to evaluate large 
mass handling situations. The mathematical model can 
also be incorporated into virtual reality training to sup- 
ply the space suit’s contribution to the virtual reality sys- 
tem’s force feedback. Physics-based models of space 
suit joint mobility will aid in designing joints for plane- 
tary space suits that have the proper combination of 
mobility and stability for best performance in walking, 
running and traversing rough terrain. An analysis com- 
bining kinematics, space suit mobility models, and 
human strength could be used to assess the lower body 
mobility of planetary space suits, relating individual 
joint mobility to terrain-traversing capabilities, such as 
the maximum height of an object that can be stepped 
over, aiding in requirement definition for advanced 
space suits. The space suit mobility database and analyt- 
ical techniques that were developed in this thesis can be 


applied to both near-term and advanced EVA planning 
and design by supplying realistic benchmarks on the 
mobility performance of the EMU space suit and 
enabling designers and planners to assess practical per- 
formance measures based on small-scale, joint level 
mobility models. 

6.6 Beyond Space Suits: Human 
Exploration Missions and Locomotion 

This research effort was enhanced by additional efforts 
to understand human performance in extreme environ- 
ments. An additional battery of experimental tests was 
run. Specifically, 

1 . The instrumented robot space suit tester was suited 
in a Shuttle EMU and commanded to perform single 
and multiple joint motions in order to collect data 
that will be used to expand the joint angle and torque 
database as well as verify the mathematical and 
physical models. 

• Suited test subjects performed representative EVA 
tasks. A Motionstar Flock of Birds motion track- 
ing system measured arm and leg joint positions 
and angles. Heart rate and surface pressure on the 
left thigh were also recorded. The tasks per- 
formed were focused around studying range of 
motion limits as well as locomotion. Tasks 
included: 

• Walking on treadmill at two different suit pres- 
sures 

• Arm reach (including a low reach test) 

• Foot locus with and with handhold 

• Task board 

With better understanding of astronaut performance we 
can look ahead to future missions and propose both 
improvements to the current space suit as well as revolu- 
tionary new designs for partial gravity locomotion space 
suits. This ongoing research is aimed at radical advances 
in our ability to understand, simulate, and predict capa- 
bilities of suited astronauts in a variety of scenarios. 

As part of the foot locus test subjects were instructed to 
make an upward spiral motion with their foot, both with 
and without a handhold. The Range of motion with the 
handhold was significantly greater than without. This 
result demonstrates the effect of balance on lower body 
mobility and suggests that an assisting device such as a 
walking stick should be utilized for human exploration 
missions (See Figure below). 
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Walking tests showed that subjects walked with longer 
steps and a lower stride frequency at lower suit pres- 
sures, giving a quantitative measure of the effect of suit 
pressure on locomotion. Stride frequencies were 
approximately 63 strides/min for 13.8 kPa (2 psi) and 70 
strides/min for 29.6 kPa (4.3 psi). 
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Funded under the NASA Institute for Advanced Con- 
cepts (NIAC), research efforts were extended to a new 
conceptual design for a space exploration system. A 
Bio-Suit System stands to revolutionize human space 
exploration by providing enhanced astronaut extravehic- 
ular activity (EVA) locomotion and life support based 
on the concept of providing a ‘second skin’ capability 
for astronaut performance 1 1 . The novel design concept 
is realized through symbiotic relationships in the areas 
of wearable technologies; information systems and evo- 
lutionary space systems design; and biomedical break- 
throughs in skin replacement and materials. By working 
at the intersection of engineering; design; medicine; and 
operations, new emergent capabilities could be 
achieved. The Bio-Suit System would provide life sup- 
port through mechanical counter-pressure where pres- 
sure is applied to the entire body through a tight-fitting 
suit with a helmet for the head. Wearable technologies 
will be embedded in the Bio-Suit layers and the outer 
layer might be recyclable. Hence, images of ‘spraying 


on’ the inner layer of the Bio-Suit System emerge, 
which offers design advantages for extreme, dusty, plan- 
etary environments. Flexible space system design meth- 
ods are slated to enable adaptation of Bio-Suit hardware 
and software elements in the context of changing mis- 
sion requirements. Reliability can be assured through 
dependence of Bio-Suit layers acting on local needs and 
conditions through self-repair at localized sites while 
preserving overall system integrity. The Bio-Suit Sys- 
tem is relevant to NASA’s strategic plan and stated 
visionary challenges in the Human Exploration and 
Development of Space, AeroSpace Technology, and 
Space Science enterprises. 

The images below depict an astronaut wearing a 
mechanical counter-pressure (MCP) suit outfitted for 
extreme Martian exploration. A ‘skin suit’ layer is com- 
plemented with climbing gear, communications, biosen- 
sors, and wearable computers for duty that requires 
high-slope traverses, repelling, and investigating craters. 
An astronaut on Mars is depicted donning the comfort- 
able elastic bio-suit layer (1). The hard torso shell (4) is 
donned next and seals with couplings at the hips and 
shoulders. The hard backpack, or portable life support 
system, (5) attaches mechanically to the hard torso shell, 
and provides gas counter pressure. Gas pressure flows 
freely into the helmet (2) and down tubes on the elastic 
bio-suit layer to the gloves and boots (3). The bio-suit 
layer is lightweight and easy to don and doff. It is cus- 
tom fitted to each astronaut using a laser scanning/elec- 
trospinlacing process (Natick Soldier Center). 
Remaining suit elements are simple, functional, inter- 
changeable and easy to maintain and repair. 
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Maintenance of a windmill within the Martian colony. 
Harvesting power from the Martian environment will be 
key to the success of Martian outposts such as this. In 
this vision we imagine windmills as one component of 
the power generation hardware. Here an individual 
scales a windmill for regular maintenance tasks. Activi- 
ties such as this push conventional space suit design to 
the limit, but this task feels natural while wearing the 
Bio-Suit. This is because of its low bulk and MCP tech- 
nologies which facilitate full range of motion as well as 
tactile feedback from the environment. Also, integrated 
into the Bio-Suit are components that allow attachment 
of such safety devices such as the Martian prussic 
(shown) or other scaling/repelling equipment. This facil- 
itates the large range of activities required of the explor- 
ers on the surface of Mars. 



not produce such fibrous material as electrospinlacing, 
but is good for creating thin elastic layers. Both melt 
blowing and electrospinlacing processes, but especially 
melt blowing, have been used in limited applications, 
commercially. Application can be made directly to the 
skin as shown here, or to advanced 3D forms generated 
by laser scanning. Wearable computers, smart gels and 
conductive materials could be embedded between poly- 
mer layers. In the illustration, melt blown gel is used to 
create a seamless MCP layer over the entire body. A 
shower-like spray device could be potentially used for 
suit self-application. 



Electrospinlacing is a process by which a multi-filament 
fiber of polymer can be projected onto a grounded sur- 
face. This is achieved by charging a suspended drop of 
polymer with tens of thousands of volts. At a character- 
istic voltage the droplet forms a Taylor cone, and a fine 
jet of polymer releases from the surface in response to 
the tensile forces generated by interaction of an applied 
electric field with the electrical charge carried by the jet. 
The projected polymer can be collected as a continuous 
web of fibers in a range of thicknesses. Application can 
be made directly to the skin as shown above, or to 
advanced 3D forms generated by laser scanning. Wear- 
able computers, smart gels and conductive materials 
could be embedded between polymer layers. In the illus- 


Advanced Technologies Melt blowing involves 
liquefying polymer and blowing it onto surfaces. It does 
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tration, electrostrictive gel is used to create a seamless 
MCP layer. 



This section introduced some related, but further reach- 
ing research efforts to enhance ongoing EVA studies. All 
illustrations in this section drawn by Cam Bresinger. 

6.7 Summary of contributions 

This research effort advances current knowledge of 
astronauts’ capabilities and limitations in space-suited 
EVA in three areas: by gathering data on the torques 
needed to bend space suit joints, developing both mathe- 
matical and physics-based models relating joint torques 
to angles, and using the joint torque-angle data to pre- 
dict and analyze the work envelope of a space suited 
astronaut. 

The space suit torque-angle database compiled in this 
thesis is more extensive in the number of joints and 
angular range than any other published space suit mobil- 
ity data set. Furthermore, the data was collected under 
realistic conditions, using joint angle trajectories that 
were generated by space-suited human subjects and 
torques measured using an instrumented robot as a sur- 
rogate for a person in the space suit. This database pro- 
vides the basis for developing and validating models of 
space suit mobility. 

The modeling work in this thesis contributes both 
numerical predictions of the torques needed to bend the 
space suit joints and insight into the physical processes 
that govern space suit joint mobility. Preisach hysteresis 
model coefficients were identified from experimental 
space suit torque-angle data and a new method was 
developed for estimating the variance of the error in the 
Preisach model’s torque predictions. The Preisach 
model was then used to generate numerical predictions 
of the torque needed to bend the space suit joint as a 
function of time. The physics-based modeling work 
compared two approximate models which describe pres- 
surized fabric cylinder bending to experimental data. 

The beam model assumes that space suit joint bending 


occurs through elastic deformations of the space suit 
fabric, while the membrane model assumes that space 
suit fabric never stretches. The experimental data agrees 
most closely with the membrane model, indicating that 
elasticity is not an important contribution to space suit 
joint bending performance and efforts to improve space 
suit joint mobility should focus on geometrical aspects 
of joint design, rather than the material properties of the 
fabric. 

The work envelope analysis demonstrates the usefulness 
of modeling space suit mobility in predicting a global 
human factors metric. The work envelope prediction 
method that this thesis developed is rapid and easily 
reconfigurable for people of different sizes and 
strengths. It generates not only boundaries on acceptable 
work sites, but also predicts the locations of desirable 
work sites. A sensitivity analysis on the work envelope 
revealed that improvements in shoulder mobility and 
upward and downward visibility would be most effec- 
tive in enlarging the space suited astronaut’s work enve- 
lope. 

Concepts for a Bio-Suit System were introduced to 
show the revolutionary potential for human space explo- 
ration. These concepts enhance astronaut extravehicular 
activity (EVA) locomotion and life support by providing 
a ‘second skin’ capability, or mechanical counter pres- 
sure for astronaut performance in extreme environ- 
ments. 
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